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The first measurements of the effect of the focal Intensity on 
the fidelity of phase conjugation in focused stimulated Brillouln 
scattering (SBS) are reported. Regardless of the strength of 
aberration, laser beams with focal intensities above 5 times the SBS 
threshold were conjugated accurately. As the focal intensity 
declined below this value, the phase conjugate fidelity declined. 

The maximum correctable aberration strength was that which reduced 
the focal intensity to near the SBS threshold. Higher laser powers 
corrected stronger aberrations. 

The pump laser was a single axial mode, tripled, Q- switched 
Neodymium:YAG, with a 6 ns pulsewidth, operating at 355 ran. The 
pulsewidth was 30 times the phonon lifetime (180 ps) of the SBS 
medium, lii^uid hexane. The I. 1 ' cm diameter beam was focused directly 
into a 2 cm diameter, 20 cm long cell using a 20 cm focal length lens 
placed 10 cm from the cell. Focal intensities ranging from below the 
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SBS threshold (22.5±4.5 GW/cm ) to 55 times the SBS threshold 
(1250±250 GW/cm^) were studied. A maximum reflectivity of 75% 
occurred at the maximum focal intensity. 

A computerized far field image analysis system directly measured 
the focal Intensity reduction due to aberrations. The focal 
intensity of the brightest far field spot controlled both the 
reflectivity and the fidelity. A pair of laminar airflow, hot wire 
array aberrators provided continuously variable, nearly sinusoidal 
phase modulation in the transverse axes. Lateral shearing 
interferometry and energy in the bucket measurements gave the Strehl 
ratio of the backscattered (reconstructed) beam, and hence the phase 
conjugate fidelity. 

A phenomenological theory for the phase conjugate fidelity of 
focused SBS was developed using resonator boundary conditions and a 
Gauss-Hermite mode decomposition. Analytical solutions for the case 
of no pump depletion were in qualitative agreement with the data. 

The first measurements of SBS gain depletion during mode beating 
within a Q-switched laser pulse are also reported. The gain 
depletion occurred only during the brightest and fastest mode beating 
spike. The remainder of the scattered pulse replicated the pump 
pulse. Gain depletion occurred when the mode beating spikes were 
approximately 2.5 times the phonon lifetime in duration. This agrees 
well with the theory of gain switching from stimulated Brillouin 
scattering to stimulated Raman scattering. 
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An undistorted laser beaa can travel long distances without much 
decrease in intensity and can be focused to a single bright spot. 

However, when a laser beaa travels through an aberrating medium, its 

phase is distorted. As the beaa continues to propagate, these phase 

distortions cause intensity distortions, turning a uniform, highly 
focusable laser beaa into a distorted, largely unfocusable beam. 

Fig. 1.1(a) is a photograph of a laser beam which has been distorted 
by passing through hot and cold flowing air, then focused onto a 
solid state video camera. It is clear that this beam has a 
nonuniform intensity structure,- and the focus consists of many spots. 

Since the energy of the original beam is shared among the spots, each 
spot has very little inte.isity. If the laser beam could be 

"reconstructed" into a single focal spot, the spot would be much more 

intense. 

One method of reconstructing distorted laser beams is optical 
phase conjugation. A "phase conjugate mirror" reverses the wavefront 
and forces the beam to travel back through the aberrator, producing & 
beam with no distortions. Such a beam is shown in Fig. 1.1(b). The 
beam is clearly more intense and sharply focused. 
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Figure 1. 






(a) 


(b) 


1 Photographs of laser far field spots. (a) - Aberrated 
spot. (b) - Phase conjugated, reconstructed spot. 
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Stimulated Brillouin scattering, Che scattering of light waves 
from sound waves, Is one of the nonlinear optical processes which can 
produce phase conjugation. Intense laser beams cause pressure 
changes (via electrostriction) in certain liquids, producing sound 
waves. These sound waves change the refractive index of the liquid, 
forming efficient scattering gratings which have the same curvature 
as the laser wavefront. The curvature of the gratings produces phase 
conjugation, or wavefront reversal of the input beam. Fig. 1.2 
depicts a distorted laser beam and the sound wave it produces. The 
sound wave produces a replica of the laser beam traveling in the 
reverse direction. 

The backward traveling wave has the same phase as the input laser 

beam, but is reversed in direction. In 1972, Zel'dovich, et. al, 

recognized that this would be characteristic of a wave which can be 

described mathematically as the the complex conjugate of the input 

wave.* - It has been shown that this process of "phase conjugation" is 

2 

equivalent to a reversal in time. The concept of time reversal is 
useful for understanding phase conjugation. A laser beam is 
aberrated as it propagates through a distorter. If this laser beam 
is time-reversed upon phase conjugation, it will automatically undo 
its aberrations as it retraverses the distorter. If the reversal is 
accurate (high fidelity), the wavefront is reconstructed to its 
original (unaberrated) quality. 
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BEAM 



SCATTERING 

MEDIUM 


Figure 1.2 Phase conjugate reflection. 







In Fig. 1.3 an Input beaa (denoted by a solid line) is aberrated as 
it passes through a distorter. It is "phase conjugated" or reversed 
in tiae by the scattering nedlua. Thus, it reflects back in the 
direction from which it caae (dashed line) and is reconstructed as it 
retraverses the aberrator. 
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Figure 1.3 Aberration and reconstruction by phase conjugation. The 


backscattered beam is automatically reconstructed. 










l.B. Background 


Spontaneous scattering of light waves by thermally excited sound 

waves was experimentally investigated by Lord Rayleigh in the 

1800's.^ In the early 1900's, Brillouin predicted that light could 

produce sound vaves by electrostriction and then scatter the sound 
4 5 

waves. ’ However, stimulated Brillouin scattering (SBS) was not 
discovered until after the invention of the laser. In the 1960's SBS 


was observed in solids, liquids, and gases 


6-10 


Even the earliest 


studies of SBS mention that the scattered beam has maximum gain in 

the backward direction** and the scattered light retraverses its input 
9 * 10 

path. 


In 1972, Zel'dovich and coworkers found that laser light which 

had passed through an etched glass distortion plates was 

reconstructed by backward SBS.* The backscattrred beam automatically 

reconstructed Itself after retraversing the distortion place. This 

unexpected result led to the concepts of phase conjugation*’* 1 and 
2 

time reversal. The fidelity, or accuracy of phase conjugation is 

generally defined as the degree of overlap of the backscattered beam 

with the input beam. The overlap can be expressed as an 
12-17 

integral. Letting F be the phase conjugate fidelity, we have 


F(z) - 


IJdxdy E L (x,y,z) E s (x,y,z)| 2 
Jdxdy |E L (x,y,z)| 2 Jdxdy|E s (x,y,z)| 1 


(1.5) 


7 








where E^ i* the Input (laser) field, Eg is the backscattered 
("Stokes") field, and z is the propagation direction. 


F(z) varies only within the Brillouin medium.^ Outside the 
medium, where F is measured, it is constant in z. If we take z-z^ 
anywhere outside the medium, 


|/dxdy E L (x,y,z Q ) E s (x.y,z Q )|» 

F - const - J- dxdy | EL (x,y,z 0 ) |a J“dxdy|E s (x,y,z 0 )| a 


( 1 . 6 ) 



The numerator is the absolute square of the Integrated overlap of the 
pump and Stokes fields, taking into account both amplitude and phase. 
The denominator is a normalization factor, which equals the product 
of the integrated pump and Stokes intensities. When F-l, the overlap 
is complete and the fidelity, or wavefront reconstruction accuracy, 
is perfect. When F-0, the overlap is nonexistent, and the wavefront 
is not reconstructed. 

Soon after the discovery and first analysis of phase conjugation, 
theoretical and experimental studies began with a goal of 
understanding the fidelity of phase conjugation. In the first 
experimental observation of phase conjugation, a long (- 1 meter), 
small cross section (- 4 mm a ), reflective tube, filled with liquid, 
was used. Because of this, the early theory and experiments 
concentrated on this "waveguide" or "lightguide" geometry. A crucial 
result of these studies is that for accurate phase conjugation, it is 
necessary to distort the beam before focusing it into the waveguide.^" 
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t 
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Theoretical studies have shown that if the bean is distorted to more 

than 10 or 15 tines the diffraction limited divergence, accurate 

phase conjugation takes place.^^^ Experimental studies have 

17 21 26-35 

confirmed this for the SBS waveguide geometry. ' 

In 1977, it was shown experimentally that accurate phase 

36 37 

conjugation occurred in focused SBS, without a lightguide. ’ In 
Ref. 37, the data indicate that the phase conjugate fidelity declined 
when the input beam divergence far exceeded the diffraction limit. 

The authors’ conclusion however, was that accurate phase conjugation 
occurs only when the divergence far exceeds the diffraction limit and 
the fidelity declines as the divergence approaches this limit. This 
is contrary to what their data show. The error is made in a 
statement analyzing the data, and in the final equation. 

This early paper was pivotal in the development of SBS fidelity 
theory, and the influence of its incorrect conclusions is still 
widespread. We note that in the common theoretical formalism, it is 
assumed that the beam must be strongly aberrated for accurate phase 
conjugation in focused SBS to occur. ^ ^ ^ ^ Many of these 
papers cite Ref. 37 as evidence for the strong aberration assumption. 

In the next several years, experimental studies of the accuracy 
of phase conjugation in focused SBS were done. None of these studies 
showed that a strong amount of aberration was required for accurate 
phase conjugation. It has been firmly established that diffraction 
limited beams are conjugated accurately by focused SBS.^^^^ 


« 
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The few studies of focused SBS with weakly aberrated beasts have shown 

accurate phase conjugation.^However, even the most recent 

theoretical works on focused SBS continue to suggest that there will 

be a minimum amount of distortion (several times diffraction limited) 

21 38 

required for accurate phase conjugation in focused SBS. ’ The 
observed maximum amounts of aberration (above which the fidelity 
declined sharply) varied widely from experiment to 
experiment.26,37,38,42 44 wide spread in the reported 

experimental data, and the different conclusions of the published 
theory and experiments have led to a lively debate about the utility 
of focused SBS for optical systems requiring accurate phase 
conjugation of weakly and strongly aberrated beams. Resolving the 
question of the minimum and maximum amounts of aberration corrected 
by focused SBS forms the crux of this study. 









This dissertation is organized as follows. The Introduction is 
followed by a description of the experimental apparatus and 
techniques. A discussion and analysis of the results follows. Next, 
a phenomenological theory of focused SBS is developed, which 
qualitatively agrees with the data. Finally, conclusions and 
recommendations for further study are presented. 

The metric system of units Is used throughout the apparatus and 
data chapters of this dissertation. Gaussian units are used in the 
theoretical chapter, because the majority of references are written 
in Gaussian units. 

There is at least one term which may be ambiguous in its usage in 
this dissertation, because of inconsistencies in common usage in the 
literature. When "pulse compression" is discussed in this 
dissertation, it generally refers to a truncation of the lowest 
intensity portions of the pulse, due to threshold effects. This is 
consistent with use in several references. In contrast, the 
compression of chirped pulses in dispersive media is not truncation. 
Rather, it is the compression of the energy of a long pulse into a 
short pulse, increasing the peak power. No such effect was observed 
in these experiments. 









2.A. Overview 


In order to measure the accuracy of phase conjugation in focused 
stimulated Brillouin scattering (SBS) over a range of aberration 
strengths, we developed an apparatus and techniques to measure the 
following parameters: 

1) The input and backscattered pulse temporal shape and spectral 
content. 

2) The input and backscattered relative intensities, which give 
the reflectivity. 

3) The input, aberrated, and reconstructed near field phases, 
giving the aberration strength and a measure of the phase 
reconstruction accuracy. 

A) The far field intensity, including the intensity profiles and 
the reduction in focal intensity caused b the aberrator. 

5) The "energy in the bucket" for the pump beam and the 


backscattered, reconstructed beam, giving the phase conjugate 
fidelity. 










These ere discussed in the following sections, but first we will 
discuss the laser, the beam train, and the SBS cell. After 
describing the diagnostics, we will discuss a phase grating optical 
aberrator, which provides easily and continuously variable 
aberrations with no optical elements to align. 







2.B. 


Recently, stimulated Brillouin scattering (SBS) has been observed 

45-49 

at ultraviolet wavelengths. In order to observe competition 

between SBS and stimulated Raman scattering, we used the tripled 

output of a Q-switched YAG laser as our pump. The laser was a Q- 

switched, line-narrowed Neodymium YAG (Quanta-Ray DCR-2A (10)). The 

tripled laser output produced up to 250 millijoules in 6 ns (FVHM) 

pulses at 355 nm. An intracavity etalon reduced the bandwidth to 
49 

less than 6 GHz (typically 1.5 GHz or better). In order to provide 

single axial mode pulses, a specially timed trigger device 

(Electronic Line Narrower, or ELN) fired the Q-switch. The ELN 

partially opened the Q-switch, allowing relaxation oscillations to 

build up. Since each relaxation oscillation underwent many passes 

through the cavity, it was significantly narrower in linewidth than a 

normally Q-switched pulse. The ELN fired the Q-switch at the peak of 

the second relaxation oscillation,^ producing single axial mode 

pulses approximately 50 percent of the time. Such techniques are 

commonly used to obtain single axial mode pulses from Q-switched 
51-54 

lasers. Using temporal and Fourier transform measurements, we 

found that the multimode pulses contained at most 5 axial modes, and 
the number of modes running varied from shot to shot. The cavity 
length L was 67 cm, giving a mode spacing of c/2L-214KHz. Thus, the 
multimode shots had less than 1.1 GHz overall bandwidth. 

The laser had an unstable resonator, high gain oscillator and a 
single pass amplifier. The laser's output energy at 1.06 microns 
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could be changed by varying Che flashlamp energy. The amount of 
chird harmonic light could be changed by rotating the doubling and 
tripling crystals away from phase matching. We used a Pellin-Broca 
prism to separate the UV light from the fundamental and second 
harmonic light. At maximum lamp energy and optimum crystal 
alignment, 250 mfllijoules of 355 nm light were produced. 

The near field and far field intensity profiles were strongly 
modulated just outside the laser. To calculate the focal Intensity 
in the SBS medium, well-defined and easily analyzed intensity profile 
was produced. The beam was apertured from the original 9 mm to 3 mm, 
and spatially filtered using a 100 cm lens and a calibrated 300 
micron diameter tungsten pinhole. The resulting beam was expanded to 
3 centimeters in diameter, and then apertured to 1.5 centimeters. 

This provided nearly constant amplitude, nearly plane wave light (the 
Strehl ratio was interferometrically measured to be greater than 
0.96). At most, 3 millijoules made it through the entire system and 
into the SBS cell. Still, this energy range allowed us to study SBS 
from threshold to 55 times threshold with accurate far field 
intensity diagnostics. 

The SBS system consisted of a 20 cm focal length, antireflection 
coated (R<0.5% per surface), plano-convex lens and a tube filled with 
HPLC grade hexane liquid. The tube, 20 cm long and 2 cm in diameter, 
was placed with the front window 10 cm from the lens. The 
reflectivity was independent of the distance from the lens to the 
cell unless the front or back window was within a centimeter of 
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focus. In Chat case, damage occurred to the window and the 
reflectivity dropped severely. The windows were "neutral solution 
treated” (Schott process) to be antireflective over a broad 
bandwidth. They had less than 2t reflectivity at 355 na when used in 
the hexane-filled cell. The windows were tilted approximately 2 
degrees from normal incidence so that specular reflections would be 
dl«card-%d. The lens and cell were arranged to collect all the 
diffracted light from the aberrator. Ve examined only the focused 
(no lightguide) geometry. 
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2.C.I. General 


A fast photodiode and transient digitizer measured the temporal 
shape of the pump and Stokes pulses. A digital Fourier transform of 
the temporal data provided the frequency spectra of both pulses. The 
following paragraphs discuss the importance of laser linevidth on SBS 
gain, and the diagnostics used to measure pump and Stokes mode 
structure, mode beating, pulse compression, and SBS reflectivity. 


There are three reasons that the gain of stimulated Brillouin 
scattering decreases as the pump laser linevidth Increases. First, 


the gain coefficient explicitly depends on the pump linevidth 


56-59 


Second, if the pump coherence length is leas than the interaction 

21 59 

length, the gain is reduced. ' Finally, mode beating can produce 

intensity fluctuations shorter than the pulse width. These intense 

56 58 59 

spikes can drive competing processes ’ ’ such as stimulated Raman 

scattering,stimulated Rayleigh scattering,stimulated 
thermal scattering,®® and optical breakdown.^® In our experiment, 
the coherence length was always longer than the interaction length. 
Mode beating did produce intensity fluctuations shorter than the 
pulse width, and the data indicates that competing processes depleted 
gain from SBS during mode beating. 
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The Brillouin gain coefficient was reduced when the laser 
operated multimode. The gain coefficient's dependence on laser 
linewidth can be written 2 ^ - ’ ^^ 


g(d*0 - 


(d*. ) 




(2.C.1) 


where C is a constant containing material parameters, is the 
laser linewidth, and Av_ is the Brillouin linewidth. For Av t «Ai/_, 
the gain is independent of the laser linewidth. When the laser 
linewidth equals the Brillouin linewidth, the gain coefficient is 
reduced by a factor of two. For the gain is inversely 

proportional to the square of the laser linewidth. 


For typical Q-svitched laser systems, the laser spectrum is 

composed of narrow bandwidth modes separated by the cavity mode 

spacing. If each mode has a narrow bandwidth, and the modes are 

spaced widely compared to the Brillouin linewidth, the SBS gain is 

59 

independent of the number of modes. The bandwidth of each mode is 

typically transform limited (approximately 27 MHz for our 6 ns 

pulses) and the mode spacing is c/2L (-214 MHz for our 67 cm cavity), 

where L is the cavity length. The Brillouin linewidth is given by 

, k 2 »> 

^B-2^r-i^“ < 2c - 2) 


where is the phonon lifetime, k is the pump wavenumber, q is the 
viscosity, and p is the density. For Hexane liquid pumped at 355 


18 










na, Che phonon lifeciae is 180 ps and Che Brillouin linevidth is 880 

MHz.^ In Chis case, each aode has a bandwidth auch less chan the 

Brillouin linewidth, but we cannot aake the assuaption that the aode 

spacing is auch greater than the Brillouin linewidth. Vith 5 laser 

aodes running, the overall laser linewidth (1.1 GHz) is slightly 

greater chan the Brillouin linewidth. The SBS gain coefficient is 

2 

reduced by a factor of l+(l.l/0.88) or 2.6. This reduction in the 
SBS gain coefficient for aultiaode shots becomes important when 
competing effects such as stimulated Raaan scattering occur. 
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2.C.2. Single Axial Mode Check 

The laser used an Electronic Line Narrower to produce single mode I 

pulses approxiaately 50 percent of the time. The Electronic Line 
Narrower, in addition to trig g ering the Q-switch, monitored the 

temporal profile of the Q-switched pulse and produced a data trigger | 

signal whenever the pulse was "smooth". The temporal pulse shape was 

also measured using an avalanche photodiode (NRC Model 877) with less 

than 200ps rise/fall time,^ and a high speed transient digitizer I 

(Tektronix Model 7912AD with 7A29 and 7B10 plugins) with less than 

68 

500 ps rise/fall time. The measurements showed that the trigger 

signal was generated whenever the Q-switched pulse had less than 5% I 

modulation. 

The spectral content of the "smooth" pulses was measured by , 

digitally Fourier transforming the temporal data. The transient 

digitizer transferred data records to an HP 320 microcomputer, which 

created a data file of intensity values as a function of time (the ^ 

sampling period was 2.56x10 ^sec). A Fast Fourier Transform program 

69 

(HP Model 98827A Waveform Analysis Software ) then analyzed this 

data file by calculating a discrete Fourier transform of the time ^ 

data and plotting the intensity as a function of frequency. Smooth 
pulses have far less power in high-frequency components than 
multimode pulses. For example, a smooth pulse (less than 5% 

» 

modulation) contains 20 times less power in the first sideband than a 
two-mode (43% modulated) pulse contains in its first sideband. The 
data trigger signal for "smooth" pulses occurred when the pulses had 

•4 
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more chan 97% of chair anergy in a single mode. This is fairly 

consisCenC with ocher reported measurements of the modulation depth 

and mode structure of Q-svitched laser pulses. Other researchers 

have predicted Chat if the ratio of the modulation depth to the 

average intensity is M, the ratio of the fundamental mode intensity 

o 

to the sideband intensity should be 16/M (see Ref. 51). For a 5% 

modulation depth, the ratio of mode intensities is predicted to be 
3 

6.4x10, and approximately 99% of the energy should be in the 
fundamental mode. 
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Since Che laser only ran single node during approximately 50 
percent of the shots, mode beating occurred during the other shots. 
The temporal and frequency replication of the pump pulse by the 
Stokes pulse was monitored during the mode beating shots. An optical 
delay was used to direct both signals into a single avalanche 
photodiode. In this way, the pump and Stokes pulses appeared on a 
single trace of the transient digitizer. The data were sent to the 
microcomputer and split into two data files. The first half of the 
trace contained the pump pulse. The second data file (the second 
half of the trace) contained the Stokes pulse. The Fourier transform 
was taken of each data file, and the Stokes spectrum was compared to 
the pump spectrum. 

For pulses well above threshold, the temporal and spectral 
content of the Stokes pulse matched the temporal and spectral content 
of the pump pulse. The mode beating data verified that the temporal 
system could respond to signals with 400 picosecond rise times and 
fall times. The Fourier transform of this data showed that the 
system could respond to signals with bandwidths as high as 1.1 GHz. 
The data acquisition system had sufficient speed and bandwidth to 
measure pulse compression and laser mode beating effects. 
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The high speed detection system measured the input and 
backscattered peak intensities, providing reflectivity data. A low 
reflection window (Schott "neutral solution" antireflection treated 
window with 1%R per surface) placed approximately one meter before 
the lens and SBS cell served as a beamsplitter (see Fig. 2.1). A 
mirror (Ml) and a detector were placed on opposite sides of the 
beamsplitter, so that the beamsplitter directed the input pulse into 
the detector, and the mirror directed the backward traveling Stokes 
pulse into the detector. The system was calibrated by replacing the 
lens and cell with a calibrated (99%R) mirror (M2) and by rotating 
the beamsplitter and the mirror Ml to give equal pump and reflected 
intensities. The mirror M2 was then removed and the lens was 
reinstalled. The SBS reflectivity was calculated from the ratio of 
the backscattered pulse peak intensity to the input pulse peak 
intensity. The lens and SBS cell window were antireflection coated, 
and passed more than 99 percent of the input and backscattered light. 
The reflectivity data were repeatable to better than 10 percent. 

It is important to note that the reflectivity data are measures 
of the intensity, not the energy. In fact, it would be possible to 
have an intensity reflectivity greater than unity if significant 
pulse compression^ ^ occurred (the energy reflectivity is always 
less than unity). However, our measurements show that pulse 
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Figure 2.1 Schematic of the time-resolved intensity diagnostics and 
optical layout. 
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compression Is only significant for pulses near threshold intensity 
(where the reflectivity is much less than unity). More will be said 
about pulse compression in the data and analysis section. 

To summarize, the time-resolved data acquisition system provided 
accurate measurements of pulse compression, mode beating, and 
reflectivity. This system provided the first direct measurements of 
subnanosecond SBS gain depletion caused by mode beating within a 
pulse. The temporal resolution was sufficient to monitor the mode 
beating and pulse compression which occurred. The system provided 
SBS reflectivity values which were repeatable to better than 10 
percent. 












Two-axis lateral shearing interferometry directly measured the 
input, aberrated, conjugated, and reconstructed wavefronts. Fig. 2.2 
is a schematic of the optical layout, emphasizing the position of the 
two shear plates. When placed before the aberrator, they could 
measure the input beam phase and the reconstructed beam phase, for 
the two beams traveling in opposite directions. When placed between 
the aberrators and the SSS cell, they measured the aberration 
strength. The interferometric measuements of the input, aberrated, 
and reconstructed beam Strehl ratios determined the phase 
reconstruction accuracy. 

Lateral shearing interferometry is discussed in Ref. 73. It is a 
self-referencing interferometric technique, combining the front and 
back surface reflections of an optical flat or wedge to give fringes 
in an overlap region which is a subaperture of the overall wavefront 
(see Fig. 2.3). The fringes generate the derivative of the wavefront 
in the direction of shear,^ so two axes of lateral shear are 
required to determine the wavefront. In order to calculate the 
wavefront from the interferograms, a VAX computer program called 
SAINT was used. The fringes for each transverse axis were digitized, 
and the data was input to the computer program. SAINT then 
calculated the wavefront, and plotted it in a pseudo 3-dimensional 
plot. The two shear axes, which were perpendicular to the 
propagation direction, were two of the axes on the plot. The optical 
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Figure 2.2 Schematic of the shearing interferometry optical layout. 
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Figure 2. 


Shear 



3 Lateral shearing interferometry. Front and back surface 
reflections interfere, giving fringes in the overlap 
region. 


28 










path difference was plotted on the third axis. SAINT also calculated 

the peak to valley and root mean square (RMS) optical path dii 'erence 

74 

and the Strehl ratio. Assuming a constant amplitude beam, if ve 

let be the RMS optical path difference (in waves) and S be the 

RMS 

Strehl ratio, then 


S - exp[ *(2*^ BM<: ) a ] 


(2.E.1) 


The wavefronts and the Strehl ratios of the input, aberrated and 
reconstructed beams gave quantitative measurements of the accuracy of 
phase reconstruction using SBS. However, these measurements are 
different from measurements of the phase conjugate fidelity F, as 
defined in chapter 1. Recall that for any Zq outside the gain 
medium, 


F - const - 


IJdxdy E^x.y.ZQ) E s (x,y,z 0 )l a 
Jdxdy |E^(x,y,z Q ) | 2 Jdxdy|E s (x,y,z 0 > | 2 


(2.E.2) 


Note that the overlap integral takes into account both the phase and 
the amplitude of the fields at Zq. In contrast, the 

interferometrically generated Strehl ratios assume constant amplitude 
beams, and only measure the effect of the phase. Although the 
interferometrically generated Strehl ratios give accurate phase 
reconstruction data, they overestimate the fidelity F because they do 
not take into account the reduction in fidelity due to amplitude 
fluctuations. 
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It was noc possible to separate the near field phase data for 
single node shots from the data for multiaode shots. Also, the time 
required to digitize the fringes (approximately 1 hour) and reduce 


the wavefront (10-15 ainutes to run, 1/2 hour to retrieve plots) 
strongly Halted the number of samples taken. Thus this technique 
was only used in the early phases of the experiment, and as a 
baseline to viois check the far field video and energy in the bucket 
measurement systems. 
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2.F. 


A computerized far field image analysis system directly measured 
the reduction in the focal intensity due to aberration. A schematic 
of the far field image analysis system is shown in Fig. 2.4. 
Beamsplitters and mirrors directed light into the lenses and camera. 
Sampling the beam between the SBS cell and the aberrator gave a 
measure of the degree of aberration, and sampling the beam after it 
was conjugated and passed back through the aberration gave a measure 
of the phase conjugate fidelity. 


The system consists of beamsplitters and mirrors (flat to A/10) 
directing light through an antireflection coated lens doublet. The 
doublet had a variable focal length, from 60 cm to over 2 m. 
Typically, the focal length was kept at 120 cm. The lenses focused 
the light directly onto the focal plane of a General Electric Model 
TN2505 Charge Injection Device (CID) array video camera. The CID 
array camera is designed to have minimum dead space between pixels. 
The far field spots produced by our laser were always three or more 
pixels wide and high, so more than 6 pixels were illuminated by every 
spot. This minimized potential problems with discrete pixels missing 
far field data. The camera was operated with an inject/inhibit 
circuit to trigger a scan for each single axial mode shot. The focal 
plane contained 242 vertical by 377 horizontal active pixels in an 
array measuring 0.66 cm by 0.88 cm, for an effective pixel size of 
27.3 nm by 23.3 pm. The camera’s response linearity (in amps per 
watt) was better than 2% for intensities ranging from 1% of 
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Figure 2.4 Schematic of the far field intensity analysis system and 


optical layout. 




saturation to 95% of saturation.^ The camera was always operated 
below the saturation intensity. Ve used neutral density filters 
placed within 10 cm of the focal plane to reduce the Intensity to 
below the camera's saturation intensity. 

A Poynting Model 505 frame grabber digitized and stored the video 
output from the camera. The frame grabber had 8-bit resolution, 
giving better than 0.4% precision. The frame grabber had memory 
capacity to store 256 horizontal by 256 vertical pixels, but since 
the camera had 242x377 pixels, there were 14 unused rows of storage, 
and 121 unstored columns of camera data. The grabber stored the 
first 256 columns of camera data, and transferred the data to a 
Hewlett Packard Model 320 microcomputer. A program called MULTI 
processed the image data and plotted the intensity profiles using a 
Hewlett Packard Model 82906A printer.^® Figs. 2.5 and 2.6 are 
typical plots of the unaberrated and the aberrated intensity 
distributions. These intensity distributions are in good agreement 
with calculated far field intensity profiles due to sinusoidal phase 
gratings (see Section 2.H.3). 

The intensity of the brightest aberrated spot (not the average 

intensity) divided by the intensity of the unaberrated spot was used 

in measurements of the reduction in focal intensity caused by 

aberrations. The brightest spot was used because it has been 

demonstrated experimentally that the intensity of the brightest spot 
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controls the SBS reflectivity. Theory predicts that the brightest 
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spot also controls the SBS phase conjugate fidelity. 


























Calibrated neutral density filters controlled the light level 
into the camera. When the beam was severely aberrated, neutral 
density filters were removed so as to increase the intensity at the 
camera. The camera's linear response gave voltage measurements that 
were proportional to the intensity. The reduction in the focal 
intensity due to aberration is 



(V .b- 8 > 

<V 85 


X N.D. 


(2.F.1) 


where I . is the aberrated intensity, 1^ is the unaberrated 
(diffraction limited) intensity, is the peak voltage level of the 
video signal for the brightest aberrated spot, is the peak voltage 
level of the unaberrated spot, B is the background voltage, and N.D. 
is the transmissivity of filters removed to observe the aberrated 
beam. We found that our laminar heated air aberrator could produce 
up to a 34-fold reduction in the focal intensity of the brightest 
aberrated spot. 


The focal intensity reduction produced by the variable aberrator 
changed by as much as 30% from shot to shot for a given 
aberrator setting, so ten shots were averaged for each aberrator 
setting. The variance of this average is /TO, or approximately a 
factor of 3 smaller than the variance of the individual measurements. 
The average value of the focal intensity reduction is therefore 
accurate to approximately 10%. These measurements were made at 
several aberration strengths, providing a wide range of focal 
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intensities. The data on the focal intensity reduction 
aberration will be presented in section 2.H. 
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2.G. 


"Energy In th 


ket" Fidelity Diagnostic 


Since energy measurements are easy to obtain and require very 
little data reduction time, they are ideal for data acquisition. A 
spatial filter and two energy meters were used to measure the 
fraction of the backscattered (Stokes) beam which was transmitted 
through the pinhole (see Fig. 2.7). For backscattered beams which 
are at least partially reconstructed (Strehl ratio greater than 0.2) 
the energy in the bucket measurements give the Strehl ratio of the 
backscattered beam. The Strehl ratio of a laser beam is defined as 
the ratio of the beam's far field intensity to the far field 
intensity of a diffraction limited beam of the same near field spot 
size and power. 


S - r 1 - (2.G.1) 

l dl 

If the pump beam is diffraction limited, the Strehl ratio of the 
Stokes beam equals the phase conjugate fidelity. Our pump beam was 
nearly diffraction limited. Its Strehl ratio was greater than 0.96. 

The energy in the bucket is a measure of the far field 
focusability of a laser beam. It is generally defined as the 
encircled energy contained within the first dark ring of an Airy 
pattern. If we define L(w) as the fraction of the energy contained 
in a circle of radius w for an Airy pattern produced by focusing a 












tripled 
yag laser 


SPATIAL 



DETECTORS 






81 

constant amplitude plane wave to the far field, we have 

L(w) - 1 - J*(kaw) - J’(kaw) (2.G.3) 

where a is the unfocused beam radius 

k is the wavevector 

J Q ^ are the Oth and 1st order Bessel functions 

This function is plotted in Fig. 2.8 (dotted line). The first point 
of inflection (at the point of intersection with the right vertical 
dashed line in Fig. 2.8) corresponds to the first dark ring in the 
Airy pattern. The energy in the bucket is the encircled energy at 
that point. For a diffraction limited beam, the energy in the bucket 
is 84 percent of the total energy. 

Compared to a diffraction limited beam, an aberrated beam has an 
increased spot size, and decreased energy in the bucket. In Fig. 

2.9, Airy patterns with varying sizes but constant integrated 
intensities are plotted. It is apparent that there is less energy in 
the bucket and lower on-axis intensity (Strehl ratio) for the 
"aberrated" spots. Fig. 2.8 shows plots of encircled energy for 
several different sized spots. We see that the bigger spots have 
less encircled energy at the bucket point, so they would transmit 
less energy through a "one Airy bucket" pinhole. 

We note that a "one Airy bucket" pinhole is not the best choice 
for Strehl ratio measurements, because the energy in the bucket is 









Figure 2.9 Normalized far field intensity profiles for several 








is not proportional to the Strehl ratio. Fig. 2.10 is a plot of the 

Strehl ratio versus energy in the bucket for a one Airy bucket 

pinhole (solid line). The data points in Fig. 2.10 were generated by 

calculating the energy in the bucket for the "aberrated" spots 

82 

represented by the curves in Fig. 2.8. It is apparent that the 
data in Fig. 2.10 do not fit a model which predicts the Strehl ratio 
co be .Linearly proportional to the energy in the bucket. 

Now, let us define the "bucket" as half the diameter of the first 
Airy ring for the unaberrated beam. A small bucket is chosen because 
it gives an encircled energy which is nearly proportional to the 
Strehl ratio (see Fig. 2.11). For a diffraction limited beam, 
approximately 60% of the energy makes it through a pinhole which is 
half an Airy bucket in diameter (left vertical dashed line in Fig. 
2.8). Aberrated beams pass proportionately less energy (see Fig. 
2.11). In order to obtain the Strehl ratio, it is necessary to 
normalize the encircled energy fraction. The normalization factor is 
1/0.6 or 1.67. A diffraction limited beam will have a normalized 
encircled energy fraction of 100%, and a Strehl ratio of 1.0. From 
this point forward, we define "energy in the bucket" as the 
normalized encircled energy fraction for a half Airy bucket pinhole. 
These measurements give the Strehl ratio of the phase conjugate beam, 
which equals the phase conjugate fidelity. 

Measurements of the phase conjugate fidelity using this energy in 
the bucket diagnostic were simple (see Fig. 2.7). Light from the 
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Energy in the Bucket 


Figure 2.10 The Strehl ratio is plotted versus the normalized energy 
in the bucket for a one Airy bucket pinhole. Note the 
nonlinear dependence. 
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Strehl Ratio 


Figure 
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Energy in the Bucket 


i 


.11 The Strehl ratio is plotted versus the normalized energy 
in the bucket for a half Airy bucket pinhole. The 
dependence is nearly linear. 
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laser vas spatially filtered by the pinhole, and beam expanded by a 
factor of 10, giving a 3 centimeter diameter central bright spot. 

This spot vas apertured to 1.5 centimeters in diameter, providing a 
nearly constant amplitude beam. The beam vas aberrated, conjugated, 
reconstructed, and focused back through the spatial filter. Because 
the beam had been apertured to half its original size, it focused to 
a spot tvice as big as the pinhole. This effectively gave a half 
Airy bucket pinhole. The beam's Strehl ratio vas simply the 
normalized energy ratio at the two detectors. The ratlometer vas 
triggered only for single axial mode pulses, so all of the energy in 
the bucket data is limited to the single axial mode case. 

The experimental normalization vas performed by directing 
apertured, diffraction limited light (verified vith shearing 
interferometry) back through the pinhole and measuring the "energy in 
the bucket". A calibrated pair of energy meters (Laser Precision 
Model Rj-7200 Energy Ratiometers) vas used to make the measurements. 
The reciprocal of the transmitted energy fraction vas the 
experimental normalization factor. This normalization factor 
depended on the reflectivity of the beamsplitters used. For 
identical beamsplitters set to identical angles, the theoretical 
normalization factor vould be 1/0.6, or 1.67. For our beamsplitters 
set to fixed angles, the experimental normalization factor vas 1.3. 
This normalization factor is accurate to approximately 10 percent. 

All the fidelity data vas multiplied by the experimental 
normalization factor. 
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The accuracy of this measurement technique can be estimated by 
examining systematic and random error contributions. From Fig. 2.11, 
it is clear that the energy in the bucket for a half Airy bucket 
pinhole is nearly linearly proportional to the Strehl ratio. The 
root mean square deviation of the points in Fig. 2.9 from the 
straight line fit is 4 percent. The experimental normalization was 
accurate to approximately 10 percent. Random errors were minimized 
by averaging 100 shots at each power and aberation setting. The 
averaged data was typically repeatable to ±4 percent. The square 
root of the sum of the squares of these errors is 11.5 percent. 

Thus, the fidelity data measured by this technique is accurate to 
better than 12 percent. 

To summarize this section, energy in the bucket measurements gave 
the Stokes beam Strehl ratio, which equals the phase conjugate 
fidelity. This technique gave rapid, accurate measurements of the 
fidelity, without extensive data reduction. 
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2.H. 


2.H.I. general 

In order Co study the fidelity of phase conjugation, it is 

desirable to use a variable aberrator. Low order Seidel aberrations 

can be produced using conventional optical machining techniques, and 

high spatial frequency random phase aberrators have been produced by 

etching glass in hydrofluoric acid. These aberrators cannot be 

varied in strength or spatial frequency without physically moving 

optical elements in the beam train. As part of a study of the effect 

of intracavity aberrations on laser beam quality, K.R. Calahan 

83 

designed and built a variable optical aberrator. Recently, ve 

analyzed the aberrator's performance,®^ comparing it to a sinusoidal 
85 86 

phase grating. ’ In this section, the aberrator design and 
performance are presented. 
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2.H.2. 


The aberrator consists of a fan which blows air through honeycomb 
channels, then over heated wires. Fig. 2.12 is a schematic drawing 
of the aberrator. A commercial 10 cm diameter "muffin fan" blows air 
into the honeycomb structures. The honeycomb material is 5 cm thick ^ 

with 0.5 cm wide hexagonal structure. «nd prvrides nearly laminar 
airflow. The air flows past nichrome wires, which are stretched in a 
grid with approximately 0.6 cm separation. The wires are heated by 1 

adjusting the current through them with a variac. The heated air 
flows perpendicular to the beam axis, producing aberrations. Two 
crossed aberrators produced aberrations in both transverse axes. The • 

beam traverses approximately 10 cm along each aberrator, accumulating 
up to 2.1 waves of optical path difference in each axis. 

A schematic of the optical layout for analyzing the aberrator is 
shown in Fig. 2.13. After spatial filtering and beam expanding, the 
collimated 355 nm light passed along both aberrators, and into the 
diagnostics. The diagnostics included lateral shearing 
interferometry and far field intensity analysis, as discussed in 
sections 2.E and 2.F. The strength and spatial frequency of 
aberration were cross checked using the near and far field diagnostic 
systems. 
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2.H.3. Aberration Data 

The near field phase data show nearly sinusoidal transverse phase 
modulation. Since our interferogram reduction software could only 
produce plots of the optical path difference (OPD) along axes at 45° 
to the aberrator axes, the number of cycles along the aberrator axes 
and their spatial frequency are reduced in the following plots by a 
factor of 72, or approximately 1.4. Bearing this in mind, 
quantitative values of the aberration strength, spatial frequency, 
and sinusoidal nature were derived from the plots. Fig. 2.14 shows 
typical plots of the OPD in one of the transverse axes for three 
different aberrator current settings. Fig. 2.15 shows sinusoidal 
curves fit to the data. The fit Is quite good: the measured OPD is 
sinusoidal to approximately 10%. By varying the current through the 
wires, the phase front was modulated by up to 2.1 waves (peak to 
valley). Table 2.1 gives the peak to valley (P-V) and root mean 
square (RMS) OPD and the interferometrically generated Strehl ratios, 
for beams passing twice through the aberrator. 
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Using the theory of sinusoidal phase gratings, the reduction in 
focal intensity due to various amounts of aberration was predicted. 
The phase grating theory predicts a roughly l/£ reduction in focal 
intensity for 2 axes of aberration, where £”0/0^ is the ratio of the 
fa* field divergence 6 to the diffraction limited far field 
divergence 0^. The reduction in the focal intensity was measured 
using the far field image analysis system described in section 2.F. 
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Figure 2.14 Plots of the measured optical path difference across the 
wavefront. 
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Table 2.1 


Aberrator OPD and Strehl Ratio Data: Double Pass 


Aberrator 1 

Voltage 

Aberrator 2 

Voltage 

OPD 

(P-V) 

OPD 

(RMS) 

Strehl 

Ratio 

0 

0 

0.15 

0.03 

0.95 

0 

10 

0.22 

0.05 

0.89 

10 

0 

0.19 

0.05 

0.91 

10 

10 

0.19 

0.04 

0.92 

0 

20 

0.64 

0.13 

0.52 

20 

0 

0.47 

0.10 

0.65 

20 

20 

0.77 

0.15 

0.40 

0 

30 

0.73 

0.18 

0.28 

30 

0 

0.56 

0.15 

0.40 

30 

30 

1.35 

0.28 

0.05 

0 

35 

1.32 

0.35 

0.007 

35 

0 

0.93 

0.22 

0.15 

35 

35 

1.35 

0.26 

0.07 
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Table 2.2 give* the OPD's, Strehl ratios, and the measured and 
phase grating-predicted values of the reduction in focal intensity 
due to a single pass along the aberrators, with selected current 
settings. The numbers in the parentheses were extrapolated from 
double pass measurements. In Fig. 2.16, the phase grating-predicted 
and measured values of the diffraction limited intensity divided by 
the aberrated intensity, and the observed divergence are plotted as 
functions of the aberrator voltage. Fig. 2.17 is a plot of the 
measured and predicted values of the focal intensity reduction of the 
brightest aberrated spot, along with a 1/f curve, plotted versus the 
measured divergence ratio (. The data is in fair agreement with the 
phase grating theory, though the intensity is greater than predicted 
for 1.5<(<10. The direct measurement of the reduction in focal 
intensity due to aberration allowed us to study, plot, and model the 
reflectivity and fidelity of focused SBS as functions of the 
aberrated focal intensity. 

In summary, the variable phase grating aberrator produced 

variable strength, nearly sinusoidal aberrations in either or both 

transverse axes. The aberrator produced up to a factor of 34 

reduction in focal intensity, in rough agreement with the sinusoidal 

8 A 8 5 

phase grating theory. ' The main drawback of the aberrator is 
that it produces relatively low spatial frequency (1.5 to 2 per cm) 
aberrations. Still, the easily and continuously variable modulation 
depth makes this aberrator useful for many optical experiments, 
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Table 2.2 


OPD, Screhl Ratio, and Predicted and Measured Focal Intensity 
Reductions: Single Pass 


Aberrator 

Voltage 

1 Aberrator 2 

Voltage 

OPD 

<P*V) 

OPD 

(RMS) 

Strehl 

Ratio 

W 1 

Pred. 

W 1 

Meas. 

«- # "dl 


0 

0 

A 

o 

(0.03) 

(0.95) 

(1) 

1 

1 

• 

10 

10 

(0.175) 

(0.04) 

(0.94) 

(1.16) 

1.05 

1.3 


20 

20 

(0.35) 

(0.09) 

(0.73) 

(3.5) 

2.25 

3.5 


30 

30 

0.72 

0.15 

.41 

11.4 

5.5 

9 

• 

40 

40 

1.65 

0.36 

.006 

43.5 

18.2 

11 


50 

50 

2.1 

0.43 

6.7x10' 

4 52.6 

28.6 

21 


60 

60 

2.1 

0.54 

10* 5 

52.6 

33.9 

24 

• 


•1 


- 
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Figure 2.17 The measured and predicted intensity reduction, and the 
reciprocal of the divergence are plotted versus the 
measured divergence. □--phase grating theory, #--data 
(1/1..), 0--1/divergence. 







including wavefront reproduction accuracy in phase conjugation. The 
sinusoidal wavefront simplifies analysis, so the aberrator is both 
easy to use and easy to analyze. 










Chapter 3. 


In this chapter, the published experimental results on the 
fidelity of phase conjugation in focused stimulated Brillouin 
scattering (SBS) are presented first. Then, results from the 
experiments which were done for this dissertation are discussed. 
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3 . 


The following effects hsve been reported In the literature 
regarding the behavior of SBS. Except as noted, these results are 
experimental. 

1) Temporal and spectral effects have been studied intensively, 

a) SBS can compress Q-switched laser pulses to less than a 


nanosecond. 


70-72 


b) If stimulated Raman scattering (SRS) has a steady state 
gain coefficient almost as large as the steady state SBS 
gain coefficient, gain switching from SBS to SRS can 


occur during short pulses. 


60-64 


87-91 

c) Mode beating can also cause gain switching from SBS 

58 59 

to SRS. ' However, there are no published reports of 
time-resolved measurements of gain switching within mode 
beating pulses. 

d) If the pump consists of a train of short pulses separated 

by less than a phonon lifetime, SBS can step up pulse by 

92 93 

pulse to a quasi-steady state value. ’ 

e) If the pump spectrum consists of narrowband spikes 

separated by more that the Brillouin linewidth 

(1/phonon lifetime), the SBS reflectivity is independent 

of the number of pump modes, and the number of Stokes 

58 59 

modes equals the number of pump modes. ’ 

f) If the pump operates broadband, the SBS reflectivity 

. 21,56,59 

declines. 
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2) If two beamlets enter the SBS medium, the reflectivity is 

79 

controlled by the focel intenelty of the brighter beam. 
Theoretical studies predict that the phase conjugate fidelity 
of the weaker beaa Is the sane as the fidelity of the 

v 80 

stronger beaa. 

3) The threshold power for focused SBS increases roughly 

37 

linearly with increasing aberration strength. 

42 

4) For low pump powers, or for pump pulses shorter than a 

38 

phonon lifetime, increasing the aberration strength 
decreases the phase conjugation fidelity, both for random 
aberrations^' 38,42,43 an< j f or astigmatism.®® None of the 
experiments measured the phase conjugate fidelity versus the 
focal Intensity. 

5) Accurate phase conjugation in focused SBS has been observed 

for beams with lS{3lOO0, where £ is the aberrated beam’s far 
field divergence divided by the diffraction limited 
divergence.28’37,38,41-44 jj on# eX p er i»e n t* showed a 

minimum divergence for accurate phase conjugation. 

6) For diffraction limited beams with focal intensities well 
above the SBS threshold, the phase conjugation fidelity 

decreases approximately 10% for a two order of magnitude 

. . 38,41-43 

increase in the pump power. 

7) Increasing the aberration strength qualitatively increases 

40 94 

the resolution of image reconstruction by focused SBS. ’ 

This may be because increasing the aberration strength 
Increases the overlap of different spatial frequencies at the 
image's focus within the SBS medium. Experiments and theory 
have shown enhanced phase conjugate fidelity when beamlets 
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in an SBS medium were aberrated in order to overlap at the 
95 96 

focus. ’ Mote that image reconstruction accuracy is 

distinct from the phase conjugation fidelity, and that Images 

can be "poorly" reconstructed even when more than 90% of the 

44 

backscattered light is conjugated. In other words, a 
nonconjugate fraction of less than 10% can cause large 
intensity modulations, giving qualitatively "poor" image 
reconstruction. 

The experiments for this dissertation were organized to resolve 
questions concerning the minimum and maximum aberration strength for 
accurate WFR in focused, steady state SBS. The experiments 
concentrate on focal Intensity effects, which have not been reported 
in the literature. The effect of the pulsewidth on the phase 
conjugation fidelity was beyond the scope of this dissertation, 
though we did investigate the influence of mode beating and pulse 
compression on the reflectivity. Image reconstruction accuracy was 
also beyond the scope of the dissertation. 
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3.B. 


3.8.1. Overview 

Our temporal and apectral data acquisition system was described 
in section 2.C. The temporal system had rise and fall times of less 
than 400 picoseconds. The Fourier transform system had a resolution 
of 100 MHz, and a nominal bandwidth of 1-1.5 GHz. We measured the 
pulse compression versus the focal Intensity, and found Increasing 
pulse compression as the intensity was reduced to near the SBS 
threshold. Ve also measured mode beating effects for pulses at 
approximately 30 times the SBS threshold Intensity. The Stokes pulse 
replicated the temporal and spectral structure of the pump pulse, 
except when the pump's mode beating intensity spikes caused competing 
processes to deplete gain from SBS. During fast, Intense mode 
beating, the Stokes pulse replicated only part of the pump pulse. 
Backward and forward stimulated Raman scattering were observed on 
some multimode shots. 
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3.B.2. Pulse Compression 


Pulse compression was studied as a function of the pump energy 
for unaberrated, single axial mode pulses. In Fig. 3.1, 
oscillographs of the pump and Stokes intensities versus time are 
presented. In the pulse shown in Fig. 3 .1(a), approximately 1.1 
millijoule entered the SBS cell, giving 25 times threshold intensity, 
approximately 0.36 mJ (1-81^) entered the cell in (b); and 
approximately 0.18 mJ (1-41^) entered the cell in (c). The SBS 
threshold occurred at 0.04510.005 mJ. The Stokes pulses were 
compressed relative to the pump pulses, with increasing compression 
as the pump energy declined. 

Table 3.1 gives the values of the pump and Stokes pulsevidths and 
compression as functions of the energy and the intensity over the SBS 
threshold. In Fig. 3.2, the pump and Stokes pulsewidths (FVHM) are 
plotted versus the input intensity divided by the threshold 
intensity. The Stokes pulse was uncompressed for a pump pulse with 
25 times threshold intensity, compressed by a factor of 0.62 for a 
pump pulse with one third that intensity, and compressed by a factor 
of 0.42 for a pump pulse with 4 times the SBS threshold intensity. 

A phenomenological model for pulse compression can be constructed 
using the intensity versus time data and the reflectivity data. 




(c) 


Figure 3.1 Pump and Stokes intensities versus time (oscillographs.) 
(a)- -1—2 51 th , (b)--I-8I th , (c)--I-4I th . 
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Table 3.1 


Pump and Stokes Pulsevidths 


Energy (mJ) 

^th 

Type 

Pulsevidth (ns) 

Compression 

1.0 

25 

Pump 

5.1 

.... 

.... 

.... 

Stokes 

5.1 

1.0 

0.34 

8 

Pump 

5.3 

.... 

.... 

.... 

Stokes 

3.3 

0.62 

0.15 

4 

Pump 

5.5 

.... 

.... 

.... 

Stokes 

2.2 

0.40 








Pump 



Stokes 




















model, Che Scokes intensity as a function of time is assumed to be 
equal to the pump intensity times the SBS reflectivity. 

In Fig. 3.3(a) the normalized pump intensity for the shot with 25 
times threshold intensity is plotted versus time (upper trace). The 
lover trace is the predicted Stokes intensity derived from the simple 
puxoe compression model. The predicted Stokes intensity has a 
maximum value of 0.62 because that is the SBS reflectivity at 25 
times threshold intensity. The predicted Stokes Intensity at each 
point equals the pump intensity times the SBS reflectivity for that 
pump intensity. The Stokes pulse in this model has FWHM-4.65 ns. 
Comparing this model to the data in Fig. 3.1(a) we see that the 
actual Stokes pulse is compressed on its leading edge, but the FWHM 
equals that of the pump. The model predicts more compression than is 
observed in this case. 

In Fig. 3.3(b) the pump and predicted Stokes intensities are 
plotted versus time for the pulse with 1-81^. The peak reflectivity 
is 0.48, and the predicted Stokes pulse has FWHM-4.15 ns for a 
compression of 0.78. The actual Stokes pulse was compressed more 
(0.62). Finally, in Fig. 3.3(c) the intensities for the pulse 
having 1—41^ are plotted. The model gives a compression of 0.62, 
while the actual compression was 0.42. 

The predicted and actual values of pulse compression are plotted 
in Fig. 3.4. The simple model is in qualitative agreement with the 
data, but it generally predicts less compression than is observed. 
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Figure 3. 


«l 



<a) 


■o 



(b) 



(c) 


•J 


3 The normalized pump and predicted Stokes pulses are 


plotted versus time. The predicted compression is 
(a)--0.92, (b)--0.78, (c)--0.62. 










Perhaps phonon buildup time effects, which were neglected here, would 
lead to additional compression because they would delay the onset of 
SBS. However, these effects would be much more difficult to model. 
They were neglected because the phonon lifetime (180 ps for hexane 
pumped at 355 nm)^ is much shorter than the pulsewidth. 
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Figure 3.4 The predicted and actual values of the pulse compression 


are plotted versus the peak pump intensity over the SBS 
threshold intensity. 






3.B.3. 


Mode beating is commonly observed in multimode laser 
87-91 

operation. Depending on the number of modes running, the 

temporal structure of a pulse may give sinusoidal or more complicated 
beating. Several distinct temporal pulse shapes were caused by mode 
beating in these experiments. The temporal pulse profiles were 
Fourier analyzed to determine their frequency structure. For up to 5 
modes running, the Stokes pulse replicates the pump with two 
exceptions: 

a) The Stokes pulse is compressed in time, and exhibits a wider 
bandwidth than the pump under single mode operation. 

b) If the mode beating of the pump leads to strong, short- 
duration intensity fluctuations at certain times in the 
pulse, the Stokes pulse has lover gain at similar times in 
its pulse. This may be due to the onset of competing 
processes. 


Pulse compression for the single mode case was discussed in 
section 3.B.2. In Fig. 3.5, the pump (left) and Stokes (right) 
pulses are plotted in time and in frequency spectrum. The pulse 
compression is obvious and the Stokes spectrum is noticeably broader. 
The Fourier transform of a narrowed temporal pulse is a broadened 
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Figure 3. 



TIME (ns) 

(a) 



FREQUENCY (GHz) 

(b) 


The pump and Stokes pulse relative intensities are 
plotted versus time and frequency. Temporal narrowing 
of the Stokes pulse leads to spectral broadening. 


5 











spectral profile. Nearly single node operation Is confirmed in the 
pump and Stokes pulses, since there is little power at the beat 
frequencies. 

The case of just two modes beating is presented in Figs. 3.6 and 
3.7. In Fig. 3.6, the beat frequency is 650 MHz, or about 3 times 
the cavity mode spacing. For the shot corresponding to Fig. 3.7, the 
beat frequency is 450 MHz, or about twice the cavity mode spacing. 
Apparently the laser can run on two modes which are not adjacent 
cavity modes. The shot in Fig. 3.6 skipped two modes, while the shot 
in Fig. 3.7 skipped one mode. 

When two adjacent modes ran, it was generally in the context of a 
total of 3 or more modes running. In Fig. 3.8, the predominant beat 
is at the cavity frequency (-> 214 MHz) but there is additional higher 
frequency modulation. 

Even for rather complicated temporal profiles, the Stokes pulse 
has nearly the same temporal and spectral structure as the pump pulse 
(see Fig. 3.9). This is in agreement with reflectivity and Fabry - 
Perot interferometry data reported in Refs. 58 and 59. Note that in 
Figs. 3.5 - 3.9, the Stokes pulse appears to have a higher intensity 
than the pump pulse. This is an artifact of the detector scheme. A 
larger fraction of the Stokes pulse entered the detector than of the 
pump pulse, so that both pulses would have approximately the same 
apparent intensity at the detector. 


76 










Figure 3. 
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Temporal and spectral structure of the pump and Stokes 
pulses for two modes beating. The modes are separated 
by 3 times the cavity mode spacing. 
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Figure 3. 








FREQUENCY (GHz) 


Temporal and spectral structure of the pump and Stokes 
pulses. Modes separated by twice the cavity mode 
spacing. 
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Figure 3.8 Temporal and spectral structure of the pump and Stokes 


pulses for adjacent modes beating. The predominant beat 
is at the cavity mode spacing, but there is other 
beating as well. 
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Multimode temporal and frequency structure. The Stokes 
pulse replicates the pump pulse. 
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For cercaln pulses, the Stokes pulse did not replicate the pump 
pulse temporally or in frequency. In Fig. 3.10 we see that the 
places in the Stokes pulse which differ from the pump pulse occur 
during the brightest mode beating spikes. These spikes are 
approximately 400 ps in duration. On the other hand, the Stokes 
replicates the pump for the rest of the pulse. This is particularly 
obvious in the upper temporal profile in Fig. 3.10. 

Although reduced SBS reflectivity due to mode beating has been 
58 59 

reported, ' to our knowledge this is the first direct observation 

97 

of partial pulse replication during mode beating. A possible 

explanation of this effect is that stimulated Raman scattering (SRS) 

56 58 59 

depletes gain from SBS in these spikes. ’ ’ In hexane pumped at 

355 nm, the steady state gain for SRS is lower than that for SBS by a 
factor of approximately 2.45 (see Refs. 56 and 62). The SRS lifetime 
is on the order of picoseconds, so it reaches steady state during the 
400 ps mode beating intensity spikes. The SBS lifetime is 180 ps, 
and the SBS gain is depleted during the 400 ps mode beating spikes. 

Ue often observed forward and backward SRS during highly modulated 
pump pulses. Second Stokes SRS shifted the wavelength of the light 
from the UV to the visible, so it was easy to see when SRS occurred. 
It seems likely that stimulated Raman scattering depleted gain from 
SBS during the most intense, fastest mode beating spikes. 
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Figure 3.10 Gain depletion in the Stokes pulse during fast, intense 


mode beating. 












3.C. 


SBS Reflectivity 

I 

The SBS reflectivity was measured as a function of the focal 
intensity for a wide range of laser power settings and aberration 
^ strengths. The reflectivity was defined to be the peak intensity of 

the backscattered pulse divided by the peak Intensity of the pump 
pulse. This quantity depended only on the focal intensity of the 
brightest aberrated spot. 


The unaberrated and aberrated relative intensities were directly 
measured, and were combined with the pulse energy and temporal data 


» 


to determine the focal intensity at each laser power and aberrator 

setting. The reflectivity was measured versus the focal intensity 

for unaberrated and aberrated pulses at various laser power settings. 

In Fig. 3.11(a), the reflectivity versus the intensity over the SBS 

threshold intensity is plotted for unaberrated pulses. The 

reflectivity exhibits a marked threshold, a monotonic increase, and a 

reduced slope as saturation became important. The maximum observed 

reflectivity was 75% at a maximum pump intensity of 55±2 times the 

SBS threshold intensity. This is in good agreement with published 

38 44 

reflectivity curves. 


The threshold intensity was precisely determined by directing the 
* backscattered light into a video camera. At threshold, backscattered 

light was observed, but Just below threshold no backscattered light 
was detected. This technique had submicrojoule responsivity. The 
> SBS threshold occurred at a focused intensity of 22.5±4.5 GW/cm’. 
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Reflectivity Reflectivity 


Figure 
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Intensity (y—) 

(b) th 


.11 The SBS reflectivity is plotted ss a function of the 

focal intensity divided by the threshold intensity. The 
intensity was changed by changing the laser power in 
(a), and by aberrating the bean in (b). 
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It is veil known and accepted, both theoretically and 
experimentally, that the product of the gain coefficient times the 
pump beam's focal Intensity times the gain length (gli) controls 
focused SBS reflectivity. ’ 37-^4,56-59 For a diffraction 

limited pump beam, the threshold value of gli is easy to calculate 
from experimental parameters. The gain coefficient g equals 0.026 
cm/MW (see Ref. 56). The peak focal intensity is the peak power 
divided by the focal spot area, times a factor of 3. The factor of 3 
is the ratio of the on-axis focal Intensity to the average focal 
intensity for a Gaussian beam. If ve let the interaction length be 
equal to a diffraction limited Rayleigh range, we have 

< 3C1 > 

where P is the peak power (Energy/FWHM for a gausslan temporal 
profile pulse), and A is the area of the focused spot. The pump 
energy at the SBS threshold was 0.045±0.005 millijoules, and the 
pulsewidth was 5.7±0.3 ns, for a power of 7.9±0.9 kilowatts. For 
A-355 nm, we have gli-17±2 at threshold. 

Now, if the gain coefficient and intensity are kept at their 
measured values, while allowing the interaction (gain) length I to be 
a free parameter, then the commonly accepted assertion that gIi-25 to 
30 at threshold^ 25,35 41 can ^e used to determine i. This makes 
the gain length equal to the following: 
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if gif-25, i-^J-1.5* (3.C. 2a) 

if gli-30 t i- 17 A “ 1,7 A (3.C.2b) 

So the gain length is 1.5 to 1.7 tines the diffraction limited 
Rayleigh range. 

Next, the beam was aberrated by various amounts at each laser 
power setting. The reflectivity (and hence the gli product) depended 
strictly on the focal intensity (see Fig. 3.11(b)--the curve is the 
same as for Fig. 3.11(a)). The data points for the reflectivity 
versus the aberrated focal intensity are in full agreement with the 
data for the reflectivity versus the unaberrated focal intensity. 
Thus, whether the beam was aberrated or the laser power was changed, 
the focal intensity controlled the reflectivity. 

Recall that the intensity of the brightest aberrated spot was 

used in these calculations (not the average intensity of the spots). 

The data show that the intensity of the brightest spot controls the 

79 80 

SBS reflectivity. Others have reported similar results. ’ 

Regardless of the amount of the beam's aberration, the 
reflectivity is controlled by the focal intensity. In Fig. 3.12, the 
reflectivity is plotted as a function of the aberration strength for 
a constant focal intensity of approximately 10 times the SBS 
threshold. We see that the reflectivity is independent of the 
aberration strength. This implies that for all amounts of 
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aberration, the gain product (gli) depends only on the focal 
intensity. Thus i, the gain length, is constant for all amounts of 
aberration. This is an important and unexpected result. 

A primary assumption of the common theoretical formalism for 
focused SBS is that the gain length is equal to the length of the 
focal constriction 13,21,37-40 anc j ot h ers 21 have found that the 
length of the focal constriction (the distance along the propagation 
direction for the aberrated beam envelope to increase in size by Jl) 
increases roughly linearly as the aberration strength increases. If 
the gain length did increase with increasing aberration strength, the 
reflectivity of an aberrated beam at a given focal intensity would 
exceed the reflectivity of an unaberrated beam at the same focal 
intensity. However, the reflectivity for unaberrated beams at a 
given focal Intensity is the same as the reflectivity of aberrated 
beams at the same focal intensity (see Fig. 3.12). Therefore, the 
gain length is the same for an aberrated beam as for an unaberrated 
beam. It is not equal to the length of focal constriction, but is 
independent of the aberration strength and equal to approximately 
1.5-1.7 times the diffraction limited Rayleigh range. 


ft 






ft 
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The longitudinal length of the interference maxima (the apeckle 

length), the transverse speckle strucure, and the amount of speckle 

crossing and twisting were measured by moving the camera along the 

focal region. These parameters were measured for beams distorted by 

the variable hot wire array aberrator (see section 2.H) and for beams 

distorted by etched glass phase plate aberrators. The phase plates 

were made by etching microscope slides in concentrated hydrofluroic 

acid for 10 to 120 seconds. They introduced strong aberrations into 

the beam (the measured increase in the far field divergence, £-0/0^, 
2 3 

was 10 to 10 ). The variable aberrator produced Ls£s24. 


The length of the speckles was nearly constant, regardless of the 
aberration strength for 1<£<1000. Fig. 3.13(a) is a plot of the 
measured values of speckle 1 ength versus the aberration strength, 
using the hot wire array aberrator, with 1<£<11. The speckle length 
is 2.0±0.8 times the diffraction limited Rayleigh range. For 
randomly etched phase plate aberrators having 100<£51000, the speckle 
length was 2.610.9 times the diffraction limited Rayleigh range (see 
Fig. 3.13(b)). Recall that the gain length was approximately 1.5 - 
1.7 times the diffraction limited Rayleigh range. This is 
approximately 0.6 times the measured values of the speckle length, 
for all aberration conditions. This leads to the conclusion that the 
gain occurs in a single speckle length in focused SBS. Note that the 
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Divergence (^) 
(b) 


Figure 3.13 The measured speckle length is plotted as a function of 
the divergence. In (a), the hot wire array aberrator 
produced 1<£<11. In (b), a series of etched glass phase 
plates produced 100<£<1000. z g c j^ i s diffraction 
limited Rayleigh range. 
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intensity decreases near the ends of each speckle, so that the gain 
length should be slightly less than che speckle length. It seems 
reasonable to assume that the gain would cut off as the Intensity 
dropped to zero at the end of a speckle. 

A primary assumption of the common theoretical formalism for 

focused SBS is that the gain length J greatly exceeds the speckle 

length f . ^^ ‘ ^ ^ In fact, it is argued that gIi B ^«l is a 

necessary condition for accurate phase conjugation in focused 
51 1R-40 

SBS. ’ ’ However our data show that t~l , and gli =25 to 

sp sp 

30, regardless of the aberration strength. This is such a novel 
result, and so directly in conflict with the published theory (that 
the existence of many speckles along the gain length is a necessary 
condition for accurate phase conjugation in focused SBS), that we 
choose to emphasize the novelty by calling it "speckle-free 
stimulated Brlllouin scattering." In contrast to the published 
theory, the gain length is not much longer than the speckle length; 
rather, it is slightly shorter than the length of the speckles. A 
large number of interference maxima along the gain length is not a 
necessary condition for accurate phase conjugation in focused SBS. 

Direct examinination of the transverse position of the speckles 
throughout the interaction length with the far field image analysis 
system, showed that the speckles do not cross or twist appreciably 
within the gain length. Even for random phase plates with £>500, the 
angle between of the speckles and the beam axisis on the order of 
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(3.D.1) 


9 

*P 


(0.5 to 2) 



where 8 is the speckle inclination angle, w is the individual 
sp sp 

speckle radius, and i is the speckle length. In other words, even 
for highly aberrated beams, the speckle crossing angle does not 
appreciably exceed the diffraction limited divergence 0^, since 


8 .i -2w /i , where w., is the diffraction limited spot size 

dl dl' Rdl sp' sp dl 

and z_is the diffraction limited Rayleigh range. 

Rdl 


The variable aberrator produced far field Intensity distributions 

which resembled the eigenmodes of a resonator. Fig. 3.14 shows 

91 

photographs of some of the modes of a stable resonator. Fig. 3.15 
shows photographs of the far field intensity distributions produced 
by the variable aberrator. Note the qualitative agreement of the 
aberrated far field intensity distributions with the Gauss-Hermlte 
modes listed below each photo. Changing the voltage in the nichrome 
wires of the aberrator changed the depth of the wavefront modulation, 
thus producing discrete far field spots. It was possible to "dial 
in" higher equivalent mode indices by Increasing the voltage. 


In the case of resonator modes, the overall beam size is nearly 
constant, regardless of the mode index. The individual spots get 
smaller as the mode index increases. However, our aberrator produces 
far field intensity patterns having individual spots which are 
constant in size w^, but having an overall size w which increases as 
the aberration strength increases. In fact, w-£w^ where { is the 
increase in the far field divergence, and w^ is the diffraction 
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Figure 3.14 Photographs of the Gauss-Hermite modes of a laser 
resonator (after Ref. 91). 
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Halted spot size. It appears that £ Is roughly equal to n+1, where 

n is the equivalent aode index. For example, the pattern eaulating *! 

the TEM^ node has an overall spot size £-n+l-4 times larger than the 

diffraction limited spot size. Note that £ equals the number of 

spots in each transverse axis, and n is the number of nulls in each 

axis. 

To summarize this section, the speckle length is Independent of 
the aberration strength, and is approximately equal to twice the 
diffraction limited Rayleigh range. The gain length is also 

independent of the aberration strength, and is approximately 0.8 J 

times the speckle length. Because the gain length is less than the 
speckle length, in contrast to current theories assuming that there 

must be many maxima and minima in the gain length, we have chosen to ^ 

■ 1 

call this "speckle-free stimulated Brlllouin scattering." The 
speckles do not cross or twist appreciably within an interaction 

length, so the equations for SBS growth should be simplified. I 

Finally, we note that using a varible phase grating aberrator, the j 

1 

focused pump intensity distribution qualitatively resembles a single I 

i 

Gauss-Hermite mode, with mode index equal to £-1, where £ is the ! 

Increase in divergence over the diffraction limited divergence. j 
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3.E.O. Overview 


tfe have found Chat Che primary influence on Che phase conjugate 
fideliCy in focused SBS is the focal intensity of the pump. 
Specifically, the intensity of the brightest aberrated spot controls 
both the SBS reflectivity and the fidelity. For intensities above 5 
times the SBS threshold, the phase conjugate fidelity, as measured by 
the energy in the bucket, is constant and equal to 0.85±0.07. As the 
focal intensity declines and approaches the SBS threshold, the phase 
conjugate fidelity declines. To first order, the effect of 
aberrating the pump is limited to its effect on the pump's focal 
intensity. As long as the aberration does not reduce the focal 
Intensity of the brightest spot to within a few times the SBS 
threshold Intensity, the phase conjugation is accurate and 
independent of the aberration strength. 


I 


I 










Shearing interferometry showed that the accuracy of phase 
reconstruction was independent of the aberration strength. Starting 
with an unaberrated intensity of approximately 20 times threshold, 
the pump beam was aberrated in either or both transverse axes. The 
phase reconstruction was exact, to within the measurement 
uncertainty, for up to 1.35 waves (peak to valley) of aberration, for 
a double pass through the aberrator. 

To find the accumulated phase introduced into an unconjugated 
beam after two passes through the aberrator, the lens and SBS cell 
were replaced with a dielectric coated mirror (R>99« at 355 run, flat 
to X/10). The wavefront phase was measured by shearing 
Interferometry, after the beam had passed tvice through the 
aberrator. Aberrated Strehl ratios as low as 0.007 were observed. 
When the lens and SBS cell were returned to the setup, the wavefront 
reconstruction was accurate to within the measurement uncertainty. 

The reconstructed Strehl ratios were all 0.92 or better, regardless 
of the strength of aberration. A Strehl ratio of 1.0 corresponds to 
a diffraction limited beam. 

In Fig. 3.16 the reconstructed Strehl ratios are plotted as a 
function of the aberration strength (input Strehl ratio). The input 
Strehl ratio was obtained from the double pass, conventional mirror 
data. In the plot, the squares are the data for one aberrator at 
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various voltage settings, the triangles represent the data for the 

other aberrator, and the circles are the data corresponding to equal 

voltage on both aberrators. The solid horizontal line at 1.0 

indicates perfect phase reconstruction. The data shows that the 

phase reconstruction in focused SBS is accurate, independent of the 

98 

aberration strength. 
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The fidelity of phase conjugation was measured using an energy in 
the bucket technique (see section 2.G. for a description of the 
technique). Strongly aberrated beams had scintillations (changes in 
intensity across the beam front), which were not removed by phase 
conjugation. Scintillations slightly decreased the conjugate beam 
near field spot size, especially for strongly aberrated beams. The 
slightly smaller near field spots were less focusable. This 
decreased the phase conjugate fidelity as measured by the energy in 
the bucket. 

In Fig. 3.17, the phase conjugate fidelity is plotted as a 
function of the aberration strength, for several pump powers. The 
fidelity monotonically declines for increasing aberration strength, 
and the rate of decline is lowest for the most Intense beam. The 
critical strength of aberration, the point where the fidelity drops 
sharply for a small increase in aberration, increases with increasing 
pump power. This appears to be consistent with the wide range of 
reported "maximum" amounts of aberration correctable by focused 
SBS.^’ 3 ^’ 3 ®’ 4 ^Higher pump powers can correct stronger 
aberrations. This data has important implications for users of SBS 
phase conjugation. The specific dependence of the fidelity on the 
focal intensity is discussed in the next section. 
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Figure 3.17 The phase conjugate fidelity is plotted versus the 
aberration strength. Before aberration, 

■~I-6I th , A~I-20I th . a-I-40I th , •~55I th . 












The effect of the focel intensity on SBS fidelity has been 
largely ignored in the published literature. Data on the fidelity 
versus the pulsevidth, Che aberration strength, and the pulse energy 
have been reported in the literature, but so far there has been no 
literature on focal Intensity effects. For this study, direct 
measurements of the reduction in the focal intensity of the brightest 
aberrated spot were made. The data show that the focal Intensity is a 
crucial parameter in the phase conjugate fidelity of focused SBS. 

The measurement technique was to set the unaberrated laser power at a 
known amount above the SBS threshold, and then increase the 
aberration strength, thereby reducing the focal intensity until SBS 
was no longer observed. Since the reduction in focal intensity for 
each aberrator setting was known, the focal intensity for each data 
point in each run was also known. 

In Fig. 3.18, the fidelity is plotted as a function of the focal 
intensity for six data runs. Within each plot, the reduction in the 
focal intensity is due to progressively increasing the aberration 
strength. Each point represents the average of 100 shots. Note that 
the fidelity is largely independent of the focal intensity, hence 
independent of the aberration strength, as long as the intensity is 
more than 5 times the SBS threshold intensity. This is particularly 
apparent in Fig. 3.18(a). The laser power was decreased for the next 
runs (b)-(f). The "critical strength" of aberration occurs at the 
point where the focal intensity is reduced to approximately five 
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Figure 3.18 The phase conjugate fidelity is plotted versus the 
aberrated focal Intensity. Before aberration, 
(a)--I-55I th , (b)--I-40I th , (c)--I-21I th , (d)--I-19I 
(e)--I-10I th , <f)-W Ith . 

10 


th’ 


3 







Claes Che SBS chreshold. In aosC cases, Che fidellcy aonoConlcally 
declines as Che aberraCion screngch increases, and as Che focal 
incenslcy decreases. 

A crucial finding of chis research is ChaC all Che daCa in Fig. 
3.18 can be flc Co a single, smooch, physically meaningful curve. In 
Fig. 3.19 all Che daca in Figs. 3.18(a)-(f) is ploCCed. The flc co 
Che curve is quite good: Che mean square devlaclon from Che curve is 
under 4 percenc. This daca was all Caken in a single afcemoon, Co 
minimize pocencial errors due Co minor alignmenc differences. The 
simple curve is physically significance Che fidelicy is constant and 
near unity (0.85 gives Che best fit) as long as the focal intensity 
is above 5 times chreshold. As the focal intensity declines to 
approach the SBS threshold, the fidelity declines sharply. 

Unaberraced and strongly aberrated data points fit the curve 
equally well. In fact, the curve in Fig. 3.19 was derived from a fit 
to the unaberrated data. The phase conjugate fidelity versus focal 
intensity data for all unaberrated shots is shown in Fig. 3.20. 

There is more spread in unaberrated intensity data than in the 
aberrated intensity data, probably because it was taken over the 
course of two weeks. Minor changes in the alignment or other factors 
could have caused the increased variance. Still, the mean square fit 
to the curve is better than 6 percent. This data is in qualitative 
agreement with the data reported in Ref. 42. The fidelity increases 
with increasing pump power when the pump is aberrated. For high 
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Figure 3.19 All fidelity versus focal intensity data from Fig. 3.18 
(a)-(f) are plotted on a single curve. 
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Figure 3.20 All fidelity versus focal intensity data for unaberrated 
pump pulses is plotted. The curve is the same as in 
Fig. 3.19. 
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pump powers, Che published decs show accurate phase conjugation, 
independent of the aberration strength. For low pump powers, 
increasing the aberration strength decreases the phase conjugate 
fidelity. 4 ^ 

To summarize this section, for both unaberrated and aberrated 

beams the intensity of the brightest focal spot controls the phase 

conjugate fidelity. The primary effect of aberration is its 

reduction in the pump beam's focal intensity. Pulses having focal 

Intensities over 5 times threshold were conjugated accurately, but 

pulses approaching the SBS threshold Intensity underwent inaccurate 

phase conjugation, regardless of their aberration strength. The 

maximum correctable aberration for focused SBS is that aberration 

which reduces the focal intensity to a few times the SBS threshold 

intensity. To our knowledge, these are the first measurements of the 

effect of the aberrated focal intensity on the fidelity of focused 
99 


SBS. 












Experiments hsve been reported that show slightly declining phase 


conjugate fidelity for unaberrated beams as the focal intensity 
increases from a few times the SBS threshold intensity to several 
hundred times the SBS threshold.^^ However, this effect is 
small: typically, a reduction of less than 10% in fidelity is 
observed for an increase in intensity of more than two orders of 
magnitude. An examiniation of the data on the phase conjugate 
fidelity as a function of the focal intensity for unaberrated beams 
reveals a similar effect. In Fig. 3.21, the fidelity versus focal 
intensity data for unaberrated (0) and slightly (£<1.3) aberrated (0) 
beams is presented. In contrast to the data presented in Fig. 3.20, 
which was taken over a two week period, this data was taken in a 
single afternoon. The line is a linear least squares fit to the 
data, indicating a 7.8% reduction in the fidelity for an Increase of 
nearly two orders of magnitude in the unaberrated focal intensity. 

No conclusive, quantitative theory has been developed to account 
for this effect, but a plausible explanation argues that the point 
where intensity exceeds threshold, hence the point where SBS occurs, 
is upstream of the focal constriction when the intensity is well 
above threshold. This means that SBS starts and phase conjugation 
occurs at other than the far field point, and not all phase 
interruptions are converted into intensity effects. The pump beam 
phase does not enter into the Stokes gain equation; only the pump 
intensity is present. Phase conjugation is accurate if the 
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Figure 3.21 The slight decrease In fidelity for increasing focal 
Intensity Is plotted for unaberrated--0 and slightly 
(£=1.3) aberrated--0 beans. The line is a linear least 
squares fit to the data. 

109 




phase fluctuations are converted to Intensity fluctuations (as 
happens at the far field point), but nay be inaccurate if there are 
many unconverted phase fluctuations (as would happen upstream of the 
focus). It would be difficult to quantify such effects, at least 
analytically. Perhaps a computer code could be developed in the 
future to take such effects into account. 

Another second order effect is a slight reduction in the fidelity 
as the aberration strength is increased, for a fixed focal intensity. 
In Fig. 3.22, the fidelity is plotted as a function of the aberration 
strength, for a fixed focal intensity of ten times the SBS threshold. 
A 10% reduction in the fidelity la observed for an order of magnitude 
increase in the aberration strength. 
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Figure 3.22 The fidelity Is plotted as a function of the aberration 
strength, for constant (1-101^) focal intensity. 
Stronger aberrations are generally conjugated with 
slightly less fidelity. 
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Using interferometrically determined Strehl ratios, for pulses 
with intensities well above the SBS threhold, the phase 
reconstruction accuracy was perfect, to within the experimental 
uncertainty, regardless of the aberration strength. Intensity 
modulation across the beam diameter was observed on the Stokes beams 
which were produced by strongly aberrated pump beams. This led to 
slightly lover values of the phase conjugate fidelity (the overlap of 
the pump and Stokes fields) than the values of the phase 
reconstruction accuracy (the pump and Stokes phase equivalence). 

Energy in the bucket measurements shoved that the aberrated focal 
intensity controlled the phase conjugate fidelity. The fidelity was 
high (0.85±0.06) and independent of the intensity unless the 
intensity was decreased to below 5 times the SBS threshold. As the 
intensity approached the SBS threshold, the fidelity dropped sharply. 
This strong dependence of the phase conjugate fidelity on the 
aberrated focal intensity may help explain the wide range of maximum 
correctable aberrations which have been reported in the literature. 
The data indicates that there will be a different maximum aberration 
for each laser power setting. Higher power pump beams can correct 
stronger aberrations using focused SBS. To the extent that it is 
possible to infer aberrated focal intensities from published 
experimental plots of the phase conjugate fidelity versus the pump 
energy and aberration strength, this data appears to be in 
qualitative agreement with published experimental results. 
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Second order effects were observed. For unaberrated beans, the 
fidelity dropped approxlnately 8* for a nearly 2 order of magnitude 
Increase in the focal Intensity. This agrees well with published 
experimental results. Also, the fidelity was consistently lower (but 
typically by less than 10 percent) for strongly aberrated beams than 
for unaberrated beams at fixed focal intensities. 
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Chapter 4. IhtgCY 


This chapter is devoted to a development of the theory of focused 
stimulated Brillouin scattering. The coupled wave equations are 
derived and solved, first using the plane wave approximation, giving 
temporal and bandwidth effects, and then using the paraxial and non- 
depleted pump approximations. A decomposition into Gauss-Hermite 
modes is used, and the equations are solved analytically for the case 
of pumping by a single Gauss-Hermite mode. This somewhat artificial 
case is qualitatively in agreement with the transverse pump Intensity 
caused by a variable phase grating aberrator. The analytic solution 
is in agreement with experimental observations that the wavefront 
reversal is accurate, and nearly Independent of the strength of 
aberration, as long as the pump intensity is well above threshold. 
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The plane wave and paraxial coupled gain equations for puatp and 

Stokes beams in stimulated Brillouin scattering are well 

known. The remarks in this section are simply a review 

of the accepted development, leading to the coupled paraxial 

equations for the pump and Stokes beams. Stimulated Brillouin 

scattering is the scattering of a pump field E^ by density- or 

pressure-driven index of refraction changes in a medium, giving a 

scattered (Stokes) field E . Specifically, the pump field scatters 

s 

from thermal phonons in the medium, giving a frequency-shifted Stokes 
field. As the number of Stokes photons increases, stimulated 
scattering can occur. The pump and Stokes waves interact to produce 
a quasi-standing wave which moves at the sound speed, instead of the 
speed of light. In steady state conditions, the Stokes wave travels 
in a direction opposite to the pump wave. The theory of Stokes wave 
buildup and its transverse intensity overlap with the pump beam is 
the subject of this chapter. 


If we Ignore absorption (the absorption coefficient is less than 
0.01 cm ^ for hexane at 355 nm) then the wave equations for a pump 
field E^ and a Stokes field E g interacting with nonlinear 
polarizations P RL(L $) are 56 


V 2 E 


L 



c 2 at 2 


a 2 p 

4s Nltl! 

c 2 at 2 


(4.A.la) 
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c 3 8t 3 


a 3 ? 

M i ku&i 

r 2 3 #- 2 


(4.A.lb) 


The flrsc Cera gives propagation and diffraction; the second gives 

temporal effects; the right hand side is the source term which 

includes bandwidth effects. The slowly varying envelope 

*4irui^ P 

approximation makes the right hand sides equal to L.S NL(L,S). 

r 2 


The nonlinear polarization is caused by pressure*induced index of 

54 

refraction changes in the medium, which are, in turn, caused by 
electrostriction. The nonlinear polarization for the pump wave can 
be written 21 * 56 ' 96 


^NL(L) * — : lfaf 2 -1- l®c' 

( ; 32* 2 cv„r B 1 + (A*/ l )V(A* b > 2 


(4.A.2a) 


where 7 is the coefficient of electrostriction, v is the velocity of 
sound, p is the density, T_ - 2xA»/_ is the Brillouin linewidth, A 1 /. 
is the laser linewidth, and Vc - + ic g is the wavevector of the 

sound wave (|k| - k * k^ - k g for backscattering). Similarly, 


P __ikl 1 _1_ |E \ 3 E 

r NL(S) 32x J cvpr B 1 + (Av l ) j /(A*/ b ) 3 


(4.A.2b) 
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The £irst tern in the right hand aide is s gain coefficient, the 

102 

second term is s Lorentzian bandwidth dependence and the third 
term is the product of fields driving the electrostriction. 


Note that r^-l/r^, the phonon lifetime. The phonon lifetime is 

. v 56 
given by 


r B 


_ e. _ 

4*1^* sin 2 (#/2) 


(4.A.3) 


where t) is the viscosity, and 4 i» the angle between the pump beam 

and the Stokes beam. The phonon lifetime becomes extremely long for 

forward scattering (4*0). Forward SBS has been predicted for pulses 

with widths more than several hundred times the backward SBS phonon 
88 

lifetime. For our 6 ns pulses, which are 30 times the backward SBS 
phonon lifetime (180 ps in hexane pumped at 355 nm),^ forward SBS 
was not observed. 





In order to simplify the notation in equations (4.A.2), we define 

_i_ 

2 “ 16*c*vpr B 1 + (A»/ L ) J /(Ai/ B ) 2 (4.A.4a) 


and 


is k7 S _ i 

2 “ 16wc J vpr_ 1 + (Av. )V<Ai/.) 2 

D JL D 


(4.A.4b) 
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Substituting the results of (4.A.4) and (4.A.2) Into (4.A.1) to 


gives 


a J E. 2i 

nip * ^ — — 

l c 2 at* 


iU fill _ HhiL |E 1 2 E 


(4.A.5a) 


an. 
sjk —S 
c 2 at 2 


■2ik c g. 


l E J ,E s 


(i.i.Sb) 


In the experiments, the 6 ns pulsewidth is 30 times the phonon 
lifetime, so ve can assume a steady state solution (exceptions will 
be discussed in section 4.C). Ve also make the paraxial 
approximation, because the fields do not vary significantly in the 

transverse direction over distances on the order of a wavelength of 

,, . CN . 16,18-26,38,39 

light. So equations (4.A.5) become 


BE 2ik T g 1 

V lh- 2Ik L^* 2 l E S l’ E L 


(4.A.6a) 
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(4.A.6b) 
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These are the coupled, paraxial, ateady state equations for 
stimulated Brillouin scattering. In the following sections we 
discuss ateady state plane wave solutions, transient solutions, and 
steady state, par-axial, undepleted pump solutions. 













4.B. 


"1 


I 


Plane Wave Equations 

I 

The solutions to the coupled equations for the pump and Stokes 
waves with no transverse derivatives are well known. They have been 
solved for the cases of a non-depleted pump,^ a depleted pump with 
no absorption, and a depleted pump with absorption.A few 
comments are in order for each of the published solutions. 


Ve shall assume a plane pump wave propagating in the +z 
direction, exciting a plane Stokes wave propagating in the -z 
direction, in a medium extending from z-0 to z-1 (see Fig. 4.1). 

The time independent equations, taking absorption into account, are 


. 'iL. , E * . a E 

8z 2 E L E S E S 2 E L 


(4.B.la) 


flE s *8c * _ 

if- fhhh + ! e s 


(4.B.lb) 


Multiplying eq. (4.B.la) by E. and adding its complex conjugate, and 


multiplying eq. (4.B.lb) by Eg and adding its complex conjugate 

. 56,100,101 

gives 


31 

3z 


"®L^L^S ” 


(4.B.2a) 
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(4.B.2b) 
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Figure 4.1 Geometry of the plane wave problem. The pump propagates 


from z-0 to z-i, but the Stokes propagates from z-i to 



I 


If we assuae no pump depletion, 1^ is s constant and we have 



‘ (S S I L * a)I S 


(4.B.3) 


This is simply exponential growth in the -z (Stokes propagation) 
direction with a gain coefficient of g c I^-a. Noting that ocO.Olcm 
for hexane pumped at 355 nm, we neglect absorption for the rest of 
the theoretical development. Thus, 


I-(0) 

Ig(J) " • x P ( «S I L i) (4B4) 

21-25 36-39 

It is widely agreed * that at threshold the gli product 

equals 25 to 30. The large amplification (10^* to 10*^) is a crucial 
factor in gain discrimination of conjugate waves from nonconjugate 
waves. Ue shall discuss the implications of intensities over 
threshold, and pump depletion below. 

If we now include pump depletion, but neglect absorption, 
equations (4.B.2) are easily integrated. We start with 



'^lVs 


(4.B.6a) 
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(4.B.6b) 
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Dividing equation (4.B.6a) by (4.B.6b) gives 


Sz g g 8z 


(4.B.7) 


so 


I L (z) 


I<;(z) - const - I (1) - 
Sc S L 


H ’ s<i> 


(4.B.8) 


If we let g-gg then I^(z)-Ig( 2 )«I^(i)-I s (i). The difference in 
intensities of the pump and Stokes waves is a constant throughout the 
gain region. Since g^/g^-w^ 8 /* 1 ^ 2 w « can think of (4.B.8) as being a 
Manley-Rowe relation,100,103 or M conservation of energy.^ The 
decrease in the pump photon flux in the +z direction equals the 
increase in the Stokes photon flux in the -z direction. 


Substituting (4.B.8) into equation (4.B.6b) we get 


' 



Sz 


- •«s I s (I s +c) 


(4.B.9) 


where C - I L * I g 


The solution is 


100,104 


I S (Z> 1 * I S (0)/I L (0) 

I s (0) " exp[(l-I s (0)/I L (0))g s I L (0)z] - I s (0)/I L (0) 


(4.B.10) 
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Unfortunately, 1^(0) is unknown. Ve can derive a transcendental 
equation for 1^(0), based on our knowledge of the noise Stokes 
intensity Ig(i) and the input pump intensity I L <0). At z-i, equation 
(4.B.10) becomes 


I s (0) exp[(l-I s (0)/l L (0))g s l L (0)i] - I s (0)/I L (0) 

I s (i) “ 1 - I s (0)/I L (0) 


I 


(4.B.11) 


and if we know the noise Stokes intensity Ig(J) then equation 
(4.B.11) is a transcendental equation in one unknown, which can be 
solved on a computer. 


If we define the reflectivity R to be I s (0)/I L (0), then equation 
(4.B.11) is simplified. 


I s (0) exp[(l-R)g s I L <0)i)-R 
I s (i) " 1-R 


(4.B.12) 


Nov, let us return to the z*dependent plane wave equatioos. 
Letting R-Ig(0)/I^(0), eq. (4.B.10) becomes 


V z > !, R _ 

I s (0) " exp[(l-R)g s I L (0)z]-R 


(4.B.13) 


We can plot equation (4.B.13) for various SBS reflectivities. At 
R-0.01, the plot of I g (z) agrees with simple exponential growth (in 
the -z direction) to better than 0.1 percent. If we assume 100% 
reflectivity (nearly 100% reflectivity has been observed with SBS in 
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For R-99%, the solution in equation (4.B.13) agrees with the 100% 
reflectivity simple form (4.B.15) to better than 5 percent. So, 
equation (4.B.13) matches the z-dependence of the simple solutions 
for conditions of low and high reflectivity. In all cases, over half 
of the total gain occurs in the first 10 percent of the gain length. 
In Fig. 4.2, we plot the z-dependent Intensity, calculated from 
(4.B.13), for several values of the reflectivity. The intensity of 
the pump monotonically declines in the +z direction, and the 
intensity of the Stokes beam is different from the pump beam by a 
constant throughout the gain region. The Stokes beam increases in 
intensity in its propagation direction (-z). 


In summary, the solutions for plane wave SBS show exponential 
growth for no depletion, hyperbolic growth for full (100%) depletion, 
and more complicated (but monotonic) growth for intermediate 
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depletion. Ve have explicitly written the transcendental equation 


(4.B.11) for the unknown constant 1^(0) in the intermediate gain 

m 

case. The gain coefficient is a lorentzian function of the laser 

linewidth divided by the Brlllouin linewidth. The gain x Intensity x 

length (gli) product is 25 to 30 at threshold. The large size of the 

gli product will be shown to give efficient, accurate discrimination I 

of conjugate waves from nonconjugate waves. 

» 

» 

•* 
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4.C. 


We have observed two temporal effects in SBS. First, for slowly 
varying pump pulses, the Stokes pulse replicates the pump pulse 
behavior, except for pulse compression. Ve have discussed temporal 
replication and pulse compression in section 3.C. Second, for 
rapidly varying pump intensities caused by mode beating, the SBS gain 
dropped and stimulated Raman scattering was observed (see section 
3.C) Competition between SBS and SRS for gain in short pulses has 
been reported elsewhere.Gain competition is analyzed below. 

Solutions to the steady state, undepleted, plane wave equations 
for SBS give 


l s - I Q exp(gli) (4.C.1) 

There are at least three sets of solutions to the transient SBS 
equations, all giving similar qualitative behavior, but slightly 
different numerical answers. In each solution the change from the 

steady state solution (4.C.1) is simply a change in the exponent. 

56 86 

For large gain, and step function pulses the gii product becomes 

gli - 2[gIi(t/r B )] 1/2 - t/r R (4.C.2) 

where t is the pulsewidth 

fg is the phonon lifetime 
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For snail gain, step function pulses, and snail laser linevldth 
conpared to the Brillouin linevldth, the gli product becomes^ 8 ’ 

gli -* gli[l * exp-(t/2fg)] (4.C.3) 


Finally, an analysis assuming constant gain throughout the 
interaction region and taking into account linewidt* elfects leads 
to 61 ’ 87 




*” ' 8li If *"ii'<'./TF 


(4.C.4) 


But if t«glifg then (4.C.4) can be reduced to 

gli - gIi(4/gIi)(t/r B ) (4.C.5) 

which is linear in t/r^. The dependence of the gain on the 
pulsewidth is plotted in Fig. 4.3, for the three models. 


We observed SBS gain reduction and visible backward and forward 

stimulated Raman scattering for certain mode beating intensity 

fluctuations with short risetimes and half widths. The phonon 

lifetime was 180 ps and the gain switching occurred for fluctuations 

having 400 ps rise and fall times, and 400 ps full widths at half 

maxima. Note that gif-25 for steady state threshold. Since g -0.026 

IS 

56 62 

cm/MW and g^-0.016 cm/MW in the steady state, ’ 
can insert measured values of the pulsewidth into equations (4.C.2-5) 


gg/gR - 2.45. We 
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3 The SBS gain reduction for short pulses is plotted versus 
the pulsewidth. The point where the gain is reduced by a 
factor of 2.5 is the Raman-Brillouin switching point. 
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and compare Che predictions of these theories to the observed gain 
switching time. 


The three curves in Fig. 4.3 bracket the observed gain switching 
time nicely. The two upper curves ought to be reduced by linewidth 
effects, making them even closer to the observed gain switching time. 
Recall that 


_ 1 

8 ‘ 1+(A* L ) 2 /(A*/ B ) a 


(4.C.6) 


But for transform limited pulses, A*/ T - l/2«rt. Also, 

la 

Av B -l/2*r B . So, 


8 " l+(r B ) 3 A J 


(4.C.7) 


This linewidth effect makes the curves approach the observed 
Brillouin-Raman gain switching condition (see Fig. 4.4). In fact, 
when the gain coefficient in eq. (4.C.2) is modified to take into 
account linewidth effects, it closely approaches the observed gain 
switching value. 
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4 The modified SBS gain reduction, taking into account 
linewidth effects, is plotted versus the pulsewidth. 
Each curve closely approaches the observed 400 ps 
switching time. 









4.D. SBS in Three Dimensions;— Phase Con1u£at lon 


The coupled paraxial equations for stimulated Brillouin 

scattering have not been solved analytically. The coupled nonlinear 

equations for second harmonic generation (an analogous process) have 

been solved approximately using a Gauss-Hermite mode decomposition 

and a perturbative approach. 1 *^ The equations for stimulated Raman 

scattering, including diffraction and depletion have been solved 

numerically, using a Gauss-Laguerre decomposition.However, 

the case of stimulated Brillouin scattering has an additional 

difficulty: the pump and Stokes waves propagate in opposite 

directions, making the boundary conditions more complex. In fact, 

109 

the coupled SBS equations are not Stuna-Liouvilie equations, so we 
have no guarantee of orthonormal eigenfunction solutions.One 
numerical approach to solving the coupled SBS equations including 
depletion and diffraction uses a Fourier decomposition of the 
transverse coordinates and an iterative "bounce" technique in the 
axial coordinate.^’ 24 ’ 2 ^ This approach may have some problems with 
stability. 24,111 


I 




If the pump is undepleted, a single equation results. Solutions 
of the Stokes gain equation have been developed for the boundary 
conditions of a lightguide (a reflective tube) when the assumption is 
made that the pump is highly aberrated. The conclusion is that for 
accurate (>0.9) wavefront reconstruction in the lightguide geometry, 
the divergence must exceed 7 to 15 times the diffraction limited 
divergence. This has been confirmed by experiments. 


I 


I 
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In a lightguide, the boundaries cause reflections, giving many 
interference maxima and minima in the gain length (see Fig. 4.5(a)). 
A. necessary condition for accurate phase conjugation in lightguide 


SBS is 


14,22 


gli g p « 1 (lightguide) (4.D.1) 


where i is the length of the interference maxima (speckles) along 
sp 

the guide. The gain along the length of the guide (gli*25) is much 
greater than the gain in any of the short interference maxima. 


For the case of focused SBS, the geometry is quite different from 
that in a lightguide. Since there are no walls to reflect the beam 
and maintain high intensity, the gain length l is shorter. Ve have 
shown in section 3.E that the gain length in focused SBS is 
independent of the strength of aberration, and is slightly less than 
the speckle length (see Figs. 4.5(b) and 4.5(c)). 

I 

In the common theoretical formalism for focused SBS, two 
assumptions are carried over from the lightguide formalism: that 

strong aberration and a gain length much longer than the speckle 

, . . , . . 13,17,21,38-40 

length are required for accurate phase conjugation. 

These assumptions are not justified. Experiments have shown that 

focused SBS produces accurate phase conjugation of diffraction 

limited beams, and with weakly aberrated beams. 

Also, we have shown that the interaction length is independent of the 

aberration strength and is slightly shorter than the speckle length. 
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Reflective Walts 




(c) 


Figure 4.5 The long gain region for SBS in a lightguide (a) is 

compared to the shorter gain region for SBS in the focused 
geometry for a diffraction limited beam (b) and for an 
aberrated beam (c). The gain length i for focused SBS is 
approximately equal to the speckle length i , which is 
! twice the diffraction limited Rayleigh range. 

f 



Not only are the fundamental assumptions incorrect, but they lead 
to incorrect conclusions. The reported solutions predict "radial 

•! 

narrowing," in which the Stokes wave's transverse envelope (cross 

section) is predicted to be approximately a factor of 2.5 narrower 

13 21 39 40 

than the pump's cross section, ’ ’ ’ nonphysical formulas for 

38 (I 

the fidelity versus the pulse duration, and inaccurate estimates of 

the minimum and maximum correctable divergence in focused 

All of the existing theories conclude that 

accurate phase conjugation is not possible in focused SBS for weakly 

aberrated beams.1^,21,37 40 j( oweveri t ^e authors of Refs. 26 and 42 

specifically mention that they observed accurate phase conjugation 

for weakly aberrated beams, contrary to theoretical predictions. * 

Similarly, in the experimental section of Ref. 40, the authors report 

that they did not observe the predicted radial narrowing. It should 

be noted that no mention is made of the aberrator's spatial frequency it 

in Refs. 26 or 42. The authors of the published theoretical papers 

generally suggest that radial narrowing and poor fidelity for weakly 

aberrated beams are not observed experimentally because pump ®. 

depletion, ignored in the theory but present in the experiments, 

would enhance the fidelity. We will show that radial narrowing and 

poor fidelity are not observed because they are outcomes of incorrect •. 

assumptions. 

There is only one paper in the literature that predicts the phase • 

conjugate fidelity as a function of the focal intensity, or as a 

13 

function of the aberration strength. After assuming a strongly 

aberrated beam (£sl0, where £ is the far field divergence over the • 
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diffraction limited far field divergence), and after assuming that 
the pump intensity profile can be approximated as a second order 
polynomial, the authors of Ref. 13 conclude that 



(4.D.2) 


In Fig. 4.6, the predicted (from Ref. 13) and measured (our data) 
fidelity is plotted as a function of the focal intensity, for (-10. 
For £>10, the predicted fidelity is even lover. There are two 
important, qualitative differences between the predictions of this 
theory and our data. First, the predicted fidelity is much lover 
than that measured, both by us and by the authors of Refs. 42-44. 
Second, the theory predicts decreasing fidelity for increasing focal 
intensity, contrary to our data and that reported in Ref. 42. The 
common theoretical formalism for focused SBS is based on Incorrect 
assumptions, and the conclusions are refuted by published 
experimental data. Therefore, we have chosen to develop a 
phenomenological theory which is not based on the incorrect 
assumptions, and which explains the data. 
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Figure 4.6 The predicted (Ref. 13) and measured phase conjugate 

fidelity is plotted versus the focal intensity over the 
SBS threshold intensity. 
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4.D.1 



We have shown in Sections 3.C and 3.D that the gain length for 
focused SBS is slightly shorter than the speckle length. This is 
reasonable because the Intensity drops near the ends of the speckles. 
The speckles did not cross or twist appreciably in a gain length. We 
have also shown that the aberrator we used produced far field 
intensity patterns that qualitatively resembled single Gauss-Hermite 
modes. To a rough approximation, the mode index n is equal to £-1, 
where £-0/0^, the divergence over the diffraction limited 
divergence. These conditions allow us to use the boundary conditions 
of laser resonators to solve the focused SBS equations. 
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We solve the equations for stimulated Brillouln scattering in a 
resonator, using a Gauss-Hermite eigenmode decomposition. We follow 



closely the analysis which was developed for stimulated Raman 
scattering, and reported in Refs. 107 and 108. Ue start with the 
undepleted, paraxial, steady-state equations: 


< 7 ! - 2l \ frh - 0 


(4.D.5a) 


and 


A ir 

< 7 ! ♦ 21k s - -1^ Ws 


where 


and 



(4.D.5b) 


85 _ l _ 

2 “ 16*c a vpr B (di/ L ) 2 /(d^ B ) 2 + 1 

- kg. Now if we solve these 
equations in two dimensions (we choose x and r) we can simplify the 
analysis. 


We note that k - k^+ kg, 


so k - |k| * k^ 
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Solutions to the homogeneous pump equation (4.D.5a) are sums of 


Gauss -Hermite functions. In x and z ve have 


E l (x,z) 


N 

- S p„ u (x,z) 


(4.D.6) 


where 


u n (x,z) “ . .1/2 , , ,-n.1/2 

((*) ' w(z)n!2 ) 


x exp[ i (kx 2 )/2R - i(n+l)tan (z/z q )] 


exp-(x 2 /w 2 (z)) H n (75x/w(z)) 


- 1 . 


(4.D.7) 


and v(z) is the beam radius, is the nth Hermite polynomial, z q is 

the Rayleigh range (*v 2 (0)/A), and R is the radius of curvature of 

the boundary: R(z) « -(l/z)(z 2 + z 2 ). We note that the pump mode 

coefficients do not vary along the z-axis. This is true for two 

reasons. First, we assume no pump depletion. Second, since the 

97 

Gauss-Hermite modes are eigenmodes of the resonator, the mode 
coefficients are constants in z. 


We can also expand the Stokes field in sums of Gauss-Hermite 
modes, making the Stokes mode expansion coefficients dependent on z. 


E s (x,z) 


A 


E 

cr—0 


S <z) u (x,z) 

O Ot 


(4.D.8) 
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Using this expansion, we eliminate transverse derivatives in the 
equation for the coefficients s^. The Stokes mode gain equation is 


‘ 11/2,8 L A A V. <W Mtp[i<-«+#)•) (*•»•«) 

n—0 m-0 0—0 


where 0 - tan*iz/z ), z q is the length of the focal constriction, 
and where 




x H (y?x/w) H<yix/w) H (7Jx/w) 
n m p 


(4.D.10) 


A variation of eq. (4.D.9) wes derived for a general mode 

expansion by Zel'dovich in the first paper to analyze phase 

conjugation.* We can think of the phase matching term exp[i(-«*+0)9] 

as a z-dependent term multiplying the gain coefficient. The real 

part decreases the gain and the imaginary part Introduces dephaslng 

between the modes s and s... We can think of the mode overlap 

a p 

integral as a four-dimensional matrix. Further assumptions 

will reduce the dimensionality of Q. Since the values of Q anm ^ are 
important to the analysis, we have calculated**^ ’ **"* and tabulated 
some Q values in Tables 4.1 and 4.2. The Q values are independent of 
ordering of the coefficients, and the largest values in each column 
or row of Table 4.1 are the "diagonal" terms The Q values are 

constants in z. 
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Values of the Overlap Integral for a-fi, n-o 


a n 

0 

1 

2 

3 

4 

0 

.398 

.199 

.150 

.125 

.109 

1 

.199 

.299 

.175 

.137 

.117 

2 

.150 

.174 

.255 

.159 

.128 

3 

.124 

.137 

.159 

.229 

.148 

4 

.109 

.117 

.128 

.148 

.210 
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T4t>l-e-6.2 








Values of the Overlap Integral Q ann ^ for a*0, n-o-0,1 


<>0001 - 0 000 

<>oin - 0 000 

<>0002 “ * -141 

<>0112 ’ 0 071 

<>0003 “ 0 000 

<>0113 “ °- 000 

w - 0 061 

<>0114 “* 092 

<>1002 - 0000 

0 m2 - 0 000 

<>1003 “ - 122 

<>1113 “ * ° 76 

w - 0 000 

<>1114 * 0 000 

<>2003 “ 0 000 

<>2113 " 0 000 

Q 2004 " * 108 

<>2114 " - 076 

<>3004 “ 0 000 

<>3114 “ 0 000 
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In Che section on boundary conditions, ve argued that to a first 
approximation, the pump field can be taken to be a single Gauss - 
Hermite mode u r throughout the interaction region (see also section 
3.E, particularly Figures 3.14 and 3.15). Increasing the aberration 
strength increased the number of spots in the transverse axis, and is 
modeled by simply increasing the mode index n of the pump. For a 
beam aberrated to have a divergence equal to £ times the diffraction 
limit, the number of spots is £, and the number of nodes n is equal 
to £-1. In this case, the pump field can be expressed as a single 
mode coefficient times the nth Gauss-Hermite mode: 

E^x.z) - p n u n (x,z) (4.D.11) 

Then equation (4.D.9) for the Stokes mode gain becomes 

» v(w>iv! ^ v <u* (4 D - 12> 

Diffraction Limited Pump 

We first solve the gain equations for a diffraction limited pump. 

For £-1, we have n-0 and the spot size w is equal to the diffraction 

limited spot size v^. The length of the focal constriction, z q , 

equals the diffraction limited Rayleigh range Zj^. The interaction 

length constrains the "resonator" boundaries to -z_ ,,<z<z_ . For 

Kal Ral 

this case, p -p_. Let us examine eq. (4.D.12), defining g' rtniI _ (l/2)g 
n u QUUp 

p o p o Q a00/r Then we have 
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(4.D.13) 


-fi¬ 

de 


s 

a 


= >MKtf V 


The mode that overlaps with the pump (the phase conjugate node) Is 
Sq. So, the equation giving the gain of the conjugate is 


-fi¬ 

de 


*0 


* *0000 *0 « P [ 1 ( * o + 0 )e] 
0—0 


(4.D.14) 


Explicitly writing out the sum, we have 

de~ s 0 “ g 0000 *0 + g 0001 *1* + *0002 *2* + < 4 - D - 15 > 

25 

We know that the total gain of s^ Is e , so this helps us deteralne 
the scale of contribution of s Q , s^, s 2 , etc., to the gain of s Q . 
Referring to our table of Q values (Table 2.2), we examine the value 
of Q 0001 and see that it is 0.000. Thus mode s^ (and in fact all odd 
numbered modes) have no contribution to the gain of s^. Qqoo 2 
-.141 (which represents loss), Qqqq 3 is 0.000, and Qqqq^ is 0.061. 

Now, even if we neglect the e term, we see that each of 

the contributions of other modes to the gain of Sq is negligible. 

(Note that we have taken exp(i(-a+£)0]-l. Using the actual value of 

exp[i(-a+0)0] can only decrease the contributions of non-phase 

matched terms.) Let us recall that the gain product gli is 

approximately 25 at SBS threshold. This product gives the intensity 

gain. The amplitude gain product is (l/2)gli, or 12.5. If we assume 

12 5 

that the contribution of mode Sq to its own amplitude gain is e , 
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we cee that the next higheet gain contributor, s^, contributes a gain 
which is related to the ratio of QqOOA^^OOOO 


contribution of s 


contribution of s, 


A _ e 12 - 5(Q 0004 /Q 0000 ) „ _-5 


12.5 


- 2.5x10 


(4.D.16) 


'id the contribution of s^ is utterly negligible in the gain of s Q . 
Ve note that for nonzero Q's their magnitude monotonically decreases 
with increasing difference between a and fi. This means that s^ and 
higher even numbered modes contribute even less to the gain of Sq 
than s^. Ve find that the only contributor to the Sq gain is itself 
(Sq), so we make a-£-0 in the gain equation (4.D.14). Gain 
selection allows vis to ignore the phase matching term exp[i(-o+/3)0] . 

The gain equation (4.D.14) for the conjugate mode s^ can now be 
written 


JL 

ae 


*o " g oooo *o 


(4.D.17) 


This gives simple exponential growth of s^. We note that the overlap 
integrals have selected only the Sq term as a contributor to its own 
gain, and non-phase matched contributions are discarded because of 
their low gain. 


Ve can use similar arguments to examine the buildup of 
nonconjugate modes. When pumped by p^, we have, as the first 
nonconjugate mode gain equation 
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*1000 "0 


e 18 ♦ ti 


1001 *1 


- g io02 s 2 * 


+ g 1003 ®3 


-216 , 
e + 


(4.D.18) 


We find Chat s^ and contribute nothing to the growth of s^, and 
the contribution of s^ (and higher odd index a^) is negligible. So, 
only contributes to its own gain. We have 

it" 8 1 “ g 1001 *1 (4.D.19) 


which again corresponds to simple exponential growth. The gain of 
the nonconjugate mode s^, compared to the gain of the conjugate, is 

to W^ooo 


g«ln of nonconjugaCe (a,) _ .12•»(Q ln01 / l Q nnn n) 
gain of conjugate (s Q ) “ e 12 ' 5 


(4.D.20) 


Thus, there is excellent gain discrimination of the conjugate mode 
from the nonconjugate. Similarly, for the next nonconjugate mode, 
Sj, we have 


and 


_ 2 _ 

ae 


- si 


2002 ®2 


(4.D.21) 


8aln of s 2 _ 

gain of s n 12.5 

0 e 


4.1x10 


(4.D.22) 
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< 


By comparing eq. (4.D.20) to eq. (4.D.22), we aee that the degree of 
discrimination against nonconjugate modes is related to the 
difference in mode indices. The larger the difference in mode 
Indices, the better the discrimination. 

Looking back over the development of the gain equations for s Q , 
s^, and when pumped by a single diffraction limited mode p^, we 
conclude that the generalised gain equation (4.D.13) can be 
simplified to 


M" *. - ®tt00a *« (4D ' 23 > 

The gain has eliminated the phase'mismatched modes. Only s 

Cl 

contributes to the gain of s^, so the equation gives simple 
exponential growth for each mode. 

We can calculate the total nonconjugate fraction of the 
backscattered field by adding up the gain fractions derived in eqs. 
(4.D.20) and (4.D.22) for the few nearest modes (the more distant 
modes contribute little to the nonconjugate, and can be neglected). 

E s 

r 0211 - Z -® (a*0) (4.D.24) 

E con ° S 0 

In order to calculate the nonconjugate fraction, we added up the gain 

112 

fractions of the first 16 nonconjugate modes. The nonconjugate 

fraction is 


149 









(4.D.25) 


r 022 - 1.9x10’ 3 + 4.1xl0* 4 + . .. - 8x10* 3 

con 

We assume that the phase conjugate fidelity equals 1-1 /I . We 

non con 

further assume that the nonconjugate Intensity fraction I /I is 

J ° non con 

equal to the nonconjugate amplitude fraction ^ non /^ con *quared. In 

*3 2 

the case of an unaberrated pump, this is equal to (2.8x10 ) , or 

7.84x10 ^. Thus, the predicted fidelity is perfect to better than 1 
part in 10 3 when the pump is diffraction limited. 

Aberrated Pump 

We have shown that for our variable aberrator, an increase in the 
divergence over the diffraction limited divergence is 

roughly equivalent to increasing the mode index (n-£-l). For 
example, if a beam is aberrated to twice the diffraction limited 
divergence, £-2 and n-1. By completely analogous arguments to eqs. 
(4.D.12-25), but for n-1, we have 

£ 

r 11211 - 1.5x10* 2 + 5.5xlO* 3 + ... - 2.3xl0’ 2 (4.D.26) 

con 

- 2 2 4 

and the nonconjugate intensity fraction is (2.3x10* ), or 5.3x10* . 
In Table 4.3, we list the nonconjugate intensity fractions for ls£<6. 
Less than 1 percent of the backscattered light intensity is 
nonconjugate, regardless of the aberration strength. 
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Nonconjugate Fraction 

as a Function of Aberration 

Strength for gIJ-25. 

(Single 

Aberration Axis) 

Divergence 

1 A 

1 /I 

non con 



•6 

1 

0 

7.84x10 ° 



* 4 

2 

1 

5.35x10 

3 

2 

9.56X10* 4 

4 

3 

1.36x10* 3 

5 

4 

1.76x10* 3 

6 

5 

2.18x10* 3 


i 


i 


i 


I 


Two Transverse Axes of Aberration 

Th« proceeding analysis examined the cese of focused SBS with a 
pump bees having intensity modulation in Just one transverse axis, 
corresponding to TEM^ nodes, where rv-{ -1. Since the equations are 
separable, ve can analyze the y-dependenca in the saae way we 
analyzed the x-dependenca. For equal aaounts of aberration in both 
transverse axes, the analysis is particularly simple. The Stokes 
mode equations take the form 

I 


L. « _ ,<* 

9B a 6 anna a 


(4.D.27) 


where 


■ ctnDtt g ■ g /% v n tt a 

* anna “ ^ ^ **n ^anna ^ 


.anna 


(4.D.28) 


and p” is the TEM mode coefficient for the pump beam 
n nn 

s° is the TEM mode coefficient for the Stokes beam 
a aa 


Thus, the two dimensional Stokes mode equations are identical in form 
to the one-dimensional equations, but the gain is proportional to the 
square of the mode overlap integral Q t U hk - This leads to improved 
discrimination of the conjugate from the nonconjugate waves for a 
given n-(-l. In other words, the maximum amount of correctable 
aberration (£ ) is higher for 2-D aberrated beams than for 1-D 

aberrated beams. 
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Up to this point, ve have assumed that the pump intensity was 
near the SBS threshold, used gIi-25, and ignored pump depletion. The 
predicted phase conjugate fidelity is nearly unity, independent of 
the aberration strength. Note that if we used gif-30, the 
nonconjugate intensity fraction would be lower and the fidelity would 
be higher. We conclude chat the phase conjugate fidelity ought to 
remain near unity for all Intensities above the SBS threshold. 


Ve can calculate the Stokes beam's amplification even during 
saturation by using our knowledge of the reflectivity as a function 
of the pump Intensity. The Stokes amplification for intensities 
above the SBS threshold is 





(4.D.27) 


where R is the measured reflectivity, is the reflectivity at 

threshold (assumed to be 0.1 percent), I L is the pump intensity, and 
I Lt k is the threshold pump intensity. So, the "effective" gain 
product gli^^ is 25+in[(R/.001)]. For intensities below 
the SBS threshold, gIJ eff *25(I L /I Lth ). 




We phenomenologically calculated the nonconjugate intensity 

fraction I /I for pump intensities above and below the SBS 
non con 

threshold, by using gll^^ instead of gIi-25 in equations (4.D.16- 
22). Even this rough phenomenological gain argument should give 
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correct qualitative predictions of the fidelity, both above and below 
the SBS threshold. Below threshold, the gli eff estimates should be 
valid. Above threshold, inaccuracies in gll^^. due to pump depletion 
simply do not matter, because the nonconjugate fraction is negligibly 
small for all intensities above the SBS threshold. 

Comparison to Experimental Data 

The values of the predicted phase conjugate fidelity are plotted 
versus the intensity over the SBS threshold in Fig. 4.7. The bold 
dashed line is the theory, and the solid line is the data. It is 
clear that the theory qualitatively agrees with the data. In the 
theory, the fidelity is near unity, regardless of the aberration 
strength, for all pump beams with focal intensities above the SBS 
threshold, and the fidelity declines for pump beams having focal 
intensities below the SBS threshold. The data shows high (*0.85) 
fidelity, regardless of the aberration strength, for all pump beams 
above 5 times the SBS threshold. The measured fidelity is slightly 
lower, and rolls off at higher intensities, compared to the predicted 
fidelity. According to this theory, the maximum correctable 
aberration for focused SBS is that aberration which reduces the focal 
intensity to below the SBS threshold. So, the maximum correctable 
aberration should be different for each laser power setting: higher 
power lasers can correct stronger aberrations. This prediction 
agrees qualitatively with the data. The measured fidelity rolled off 
as the aberration reduced the focal intensity to approximately 5 
times the SBS threshold. 
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Figure 4.7 The phase conjugate fidelity Is plotted as a function of 
intensity over the SBS threshold. The solid line is 
a fit to our data. The dashed line is our theoretical 
prediction. 
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There are several possible reasons why the measured phase 
conjugate fidelity is slightly worse than our theory predicts. The 
many assumptions that were used in our theory to make it analytically 
tractable may have led to solutions which give predictions of better 
fidelity than is observed. In particular, the assumption of the pump 
being a single Gauss-Hermite mode was a very rough approximation. 

The pump intensity profile gave an appropriate number of spots in 
approximately the correct locations, but the intensity of the spots 
decreased, rather than increased (as do the Gauss-Hermite modes) near 
the edges of the pattern. Also, some of the speckles ended or began 
in the gain region. Both these observations show that using a single 
Gauss-Hermite mode for the pump is a crude approximation. 

It is clear that more work needs to be done to make the theory 
more complete and accurate, especially for near-threshold 
intensities. A full numerical investigation of the propagation and 
gain of the pump and Stokes beams was beyond the scope of this 
dissertation, but such a project might be appropriate for a 
theoretical dissertalon. Perhaps the work in Refs. 107 or 108 could 
be extended to include a Gauss-Hermite mode decomposition. 

To summarize this chapter, a theory for the phase conjugate 
fidelity of focused SBS has been presented. A Gauss-Hermite mode 
decomposition was used in the undepleted pump limit, and the pump was 
assumed to be a single Gauss-Hermite mode. In contrast to existing 
theories, which predict inaccurate phase conjugation for weakly 
aberrated beams, this theory predicts accurate phase conjugation, 
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independent of the aberration strength, for all intensities above the 
SBS threshold. This is in qualitative agreement with the data. 



Chapter 5. 


Ue have experimentally measured the effect of the aberrated focal 
intensity on the stimulated Brillouin scattering (SBS) phase 
conjugate fidelity. We directly measured the reduction in focal 
intensity of the brightest aberrated spot and used this intensity 
(not the average focal intensity) in our analysis of the data. 
Regardless of the aberration strength, if the focal intensity was 
greater than 5 mes the SBS threshold intensity, the fidelity was 
high (greater than 0.8). As the focal Intensity was reduced (either 
by aberrating the beam or by turning down the laser power) so that 
the focal intensity approached the SBS threshold intensity, the 
fidelity declined. To our knowledge, this is the first direct 
measurement of the effect of the aberrated focal intensity on SBS 
phase conjugate fidelity. 

These measurements appear to be in qualitative agreement with 
published experimental results. In some published experimental 
results, the fidelity improves as the laser power increases, for a 
given aberration strength. Other data shows that the fidelity 
improves as the aberration strength decreases, for a given laser 
power. Our data may explain the wide variation in the reported 
values of the maximum correctable aberration strength for focused 
SBS. We obtained a different value of the "maximum" correctable 
aberration strength for each laser power setting. 


158 











The measurements refute a primary assumption of the common 
theoretical formalism: that strong aberration is a necessary 
condition for accurate phase conjugation in focused SBS. Ve found 
that the fidelity decreases for increasing aberration strength, 
because the aberration reduced the focal intensity. Our reflectivity 
measurements showed that the gain length was independent of the 
aberration strength, and was slightly less than the length of the 
interference maxima (speckles) at the focus. This data refutes 
another assumption of the common theoretical formalism: that many 
speckles must occur along the gain length for accurate phase 
conjugation in focused SBS. Also, our data refute the conclusion in 
the published theory that the phase conjugate fidelity would be 
sharply reduced due to a significant radial narrowing of the 
backscattered beam compared to the pump beam. The radius of the 
backscattered beam was approximately equal to the radius of the pump 
beam. 

Since our data refuted the assumptions and the conclusions of the 
common theoretical formalism, we developed a phenomenological theory 
which did not use those assumptions. Ve assumed that the equations 
for focused, undepleted SBS could be solved using the boundary 
conditions of a laser resonator. By using a Gauss-Hermite mode 
decomposition, and by assuming that the pump could be described by a 
single Gauss-Hermite mode, we were able to derive simple analytic 
expressions for the gain of the conjugate and the nonconjugate modes. 
We calculated predictions of the nonconjugate intensity fraction of 
the backscattered light for several aberration strengths and spatial 











frequencies. The theory predicted accurate phase conjugation, 
independent of the aberration strength, for all intensities over the 
SBS threshold. The theory qualititively agrees with our data and 
other published experimental results. 

Ve also measured the accuracy of temporal and spectral 
replication of the pump pulse by the backscattered (Stokes) pulse. 
Using time-resolved intensity diagnostics and Fourier analysis, we 
found accurate temporal and spectral replication, except during pulse 
compression and intense mode beating. Ve developed a 
phenomenological model for the pulse compression based on observed 
values of the reflectivity as a function of Intensity. The model was 
in qualitative agreement with the pulse compression data. During 
mode beating, the Stokes pulse replicated the punrp pulse, except for 
the most Intense spike. During that spike, SBS gain depletion 
occurred. To our knowledge, this is the first time-resolved 
measurement of partial pulse replication by SBS. Ve assumed that 
stimulated Raman scattering depleted the SBS gain during mode beating 
spikes. Ve found that the pulsewidths for which we observed SBS gain 
depletion wer~ as predicted in the published theory. 

Ve developed two new techniques and devices for this 
dissertation. Ve used a variable phase grat.'->g aberrator to distort 
the pump beam, and we directly measured the focal intensity reduction 
due to the aberrator, using a computerized far field video image 
analysis system. These techniques and devices may be useful for 












other researchers studying phase conjugation and imaging accuracy in 
a variety of nonlinear optical processes. 

Perhaps the most intriguing unresolved question has to do with 
image reconstruction accuracy. A quantitative study of the spatial 
frequencies of images reproduced by focused SBS under a vide range of 
aberration conditions would make an excellent follow-on project. 
Another fruitful area for further study would be a computer analysis 
of focused SBS, in the undepleted pump regime, when the pump beam 
consists of a sum of Gauss-Hermite modes (Instead of a single mode). 
Perhaps an existing code which has been developed to study stimulated 
Raman scattering using a Gauss-Laguerre mode decomposition could be 
modified to analyze non-rotationally symmetric intensity 
distributions, if the Gauss-Hermite decomposition were used. 

Finally, an experimental study of the phase conjugate fidelity of 
focused SBS as a function of the aberration strength for irregular or 
turbulent aberrations would be an attractive point of comparison with 
the regular, smooth aberrations studied in this dissertation. The 
simplicity and applicability of SBS for many uses of phase 
conjugation ensures that much more research will be done on this 
fascinating physical effect. 
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